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ABSTRACT 

Using the United Kingdom Infrared Telescope on Mauna Kea, we have carried out 
a new near-infrared J, H, K monitoring survey of almost a square degree of the star¬ 
forming Orion Nebula Cluster with observations on 120 nights over three observing 
seasons, spanning a total of 894 days. We monitored ~ 15,000 stars down to J Ri 20 
using the WFCAM instrument, and have extracted 1203 significantly variable stars 
from our data. By studying variability in young stellar objects (YSOs) in the H — K,K 
color-magnitude diagram, we are able to distinguish between physical mechanisms of 
variability. Many variables show color behavior indicating either dust-extinction or 
disk/accretion activity, but we find that when monitored for longer periods of time, 
a number of stars shift between these two variability mechanisms. Further, we show 
that the intrinsic timescale of disk/accretion variability in young stars is longer than 
that of dust-extinction variability. We confirm that variability amplitude is statistically 
correlated with evolutionary class in all bands and colors. 

Our investigations of these 1203 variables have revealed 73 periodic AA Tau type 
variables, many large-amplitude and long-period {P > 15 day) YSOs, including three 
stars showing widely-spaced periodic brightening events consistent with circumbinary 
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disk activity, and four new eclipsing binaries. These phenomena and others indicate the 
activity of long-term disk/accretion variability processes taking place in young stars. 
We have made the light curves and associated data for these 1203 variables available 
online. 

Subject headings: stars: formation - stars: pre-main sequence - 


1. Introduction 


The Orion Nebula Cluster (ONC) is the richest young cluster in the solar vicinity (for detailed 
Muench et al.|[2008 and O’Dell et al.||200^ ). At a distance of only ~414±7 pc (Menten 


reviews, see 


et al. 2007|, it is possible to study the cluster in great detail. The classic study of stars towards 


the ONC is due to Parenago (1954), followed by the work of Herbig &: Terndrup (1986), Jones & 


Walker (1988), and Hillenbrand (|1997), and most recently by catalogues of stars based on optical 


and near-infrared HST data ( Da Rio et al.|[2009 2010, 2012; Robberto et al.|[2013 ). Major progress 


has been made in identifying the population of young stars at various evolutionary stages in the 
ONC. Large numbers of Ha emission line stars have been found centered on the Orion Trapezium 


(e.g-, 

Haro et al. 1953 

Parsamian Sz Chavira 1982, Wiramihardja et al. 

1991 

Szegedi-Elek et al. 

2013 

Pettersson et al.|2014). Membership from near- and mid-infrared studies has been determined 

by, e.g.. 

Hillenbrand Sz Carpenter (2000), Lada et al. (2000) 

Muench et a 

1. (2002) 

Robberto et al. 

(2005 

), 

Andersen et al. 

(2011), Megeath et al. (2012D, and 

Billot et al. 

(2012 

). Similarly, X-ray 


2008, Pillitteri et al.|[2013 ). 


The first variable star recognized in the ONC, now known as AF Ori, was found around 1848 


by W.C. Bond (Herbig 1982). During the 1920s and 1930s, attention was drawn to numerous ir¬ 


regular variables in the region of the Orion Nebula, and surveys to discover and characterize these 


stars were carried out (e.g., Pickering &: Leavitt 1904 Ludendorff 1928 [Brun 1935, Hoffmeister 


1949, Parenago 1954, Rosino 1956). A variety of classification schemes based on light curve be¬ 


havior were devised (for further references, see the monographs by Glasby 1974 and [Hoffmeister, 


Richter, &: Wenzel 1984), but none provided an important key to understand the nature of these 


stars. The spectroscopic identification of T Tauri stars as a class (Joy 1945, Herbig 1962) and 


the recognition of T Tauri associations (Ambartsumian 1949) eventually provided the basis for a 


physical understanding of the large populations of variable stars associated with dark clouds. 

The photometric behavior of young stars in associations is best described as irregular, with 
an enormous range in amplitude, time scale and frequency of fluctuations. This is discussed in 
more detail in subsequent sections. It gradually became evident that the variability of young stars 
is difficult to characterize, and in particular that the type of variability, e.g., the RW Aur stars, 
does not define homogeneous groups of stars (e.g., Herbig||1960 ). Spectroscopy, on the other hand, 
emerged as a powerful tool to probe the nature of young stars, reducing interest in photometric 
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studies. In recent years, wide-field optical and infrared array detectors have become much more 
common, and as a result large-scale photometric surveys have again attracted interest. 

Optical surveys have mainly focused on periodic variables, due to cool or hot spots on the 


stars (e.g., Stassun et al. 1999 Herbst et al. 2000, 2002; Rodn'guez-Ledesma et al. 2009). A few 


studies have also focused on the irregularly variable stars (e.g., Parihar et al. 2009), and lately 
infrared studies, including those performed with Spitzer and Herschel, have produced a wealth of 


results (Carpenter et al. 2001, Morales-Calderon et al. 2011, Billot et al. 2012). Carpenter et al 


(2001) studied a 0?8 x 6° strip in the Orion A molecular cloud centered on the ONC and complete to 
Kg = 14.8. This ~ 50-day study showed variability arising from starspots, dust extinction, accretion 
variability, and changes of the inner disk. However, no previous infrared monitoring campaign in 
the ONC has had more than two months of continuous coverage, which limits our understanding of 
long-term variability, especially if young stars change their short-term behavior over long timescales. 
Strong suggestions of such long-term behavior were seen in the JHK monitoring study of a small 


number of YSOs in the Cygnus OB7 star-forming region presented by Rice et al. (2012) and Wolk 


et al. (2013). 


We have carried out monitoring observations of the Orion Nebula Cluster in the J, H, K near- 
infrared bands on 120 nights spanning 894 days from late 2006 to mid-2009. These observations, 
obtained using the WFCAM instrument on the United Kingdom Infrared Telescope (UKIRT), have 
allowed us to identify 1203 variable stars within a 53 x 53 arcmin region centered on the Trapezium 


stars. The evolutionary class of nearly all of these 1203 variables has been determined by Megeath 


et al. (2012), allowing us to derive NIR variability properties as a function of class. As established in 


previous JHK surveys of young stars (cf. [Carpenter et ar]|2001 Rice et aT]|2012 Wolk et ar]|2013 ), 
different variability mechanisms impart different color signatures m. JHK data: starspot-induced 
variability is nearly colorless, variable dust extinction makes stars redder as they get fainter, and 
accretion variability or other changes in inner disk geometry cause stars to get redder as they get 
brighter. We identify these different motions in H — K,K color-magnitude space and measure the 
fraction of stars showing either dust-reddening or disk/accretion variability, noting that the fraction 
of disk/accretion variables increases with longer monitoring campaigns. 

The UKIRT observing campaign is described in Section along with the data reduction. Our 
strategy for selection of variables and other methodology is described in Section to select vari¬ 


able stars we use the Stetson index (Stetson 1996); we identify periods through the Lomb-Scargle 


periodogram ( Press fc: Rybicki|l989 ) and the fast chi-squared algorithm (fy^, Palmer|[2009 ); and we 
cross-match many of our targets to catalogs in the literature in order to obtain further information, 
such as the Spitzer-devived evolutionary classes from Megeath et al. (2012). In Section]^ we show 
that variability mechanisms can be constrained via the color-color and color-magnitude diagrams; 
we derive the fraction of stars that show blueing behavior (associated with accretion/disk activ¬ 
ity) versus reddening behavior (associated with dust occultation), and find that this ratio evolves 
significantly with the monitoring duration; many individual systems were seen to switch between 
color variability modes during this campaign. In Section we show that different color variability 
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modes have different timescales, with blueing (accretion/disk activity) variability showing longer 
typical timescales. We report a number of exceptional variables in Section including a newly- 
identified class of stars with widely-separated periodic brightening events, “dipper” objects (cf. 
Morales-Calderon et al. 2011[ [Stauffer et al.|[^014) and periodic AA Tau variables, and other stars 


with long IR periods; Section identifies four previously unreported eclipsing binary systems. We 
connect our results with other recent infrared monitoring surveys in the concluding commentary in 
Section [H 


2. OBSERVATIONS AND DATA 


2.1. Observations: The UKIRT Campaign 


Our data comprise J, H, K observations of a 0.9 x 0.9 degree field towards the Orion Nebula 
Cluster that were taken with the Wide Field Camera (WFCAM) instrument on the 3.8 m United 
Kingdom Infrared Telescope (UKIRT), which sits atop Mauna Kea, Hawaii at 4194 m elevation. 
Observations have been taken on 120 nights between October 2006 and April 2009 over three 
observing seasons: October 26, 2006 -- April 21, 2007; November 21, 2007 - December 2, 2007; 
and November 8, 2008 - April 7, 2009. This time coverage is longer than any previous infrared 
monitoring campaign in the ONC; see Figure for a time coverage comparison. 


The WFCAM instrument, described in jCasali et al. (2007), consists of four Rockwell Hawaii-II 
PACE arrays of 2048 x 2048 pixels. The detectors are widely spaced such that four pointings are 
required to image a contiguous square field on the sky (Figure]^; this 2x2 observing strategy 
ultimately defines the 53' x 53' field of view of our survey. These ONC observations use the same 
instrument and configuration outlined in Rice et al. (2012, see that paper for a more detailed 


description of this observing setup). WFCAM was first used to study YSO variability by Alves de 


Oliveira & Casali (2008). 


The WFCAM observations are centered on the coordinates a 2000 = 05:35:21.0, <52000 = —05:24:35, 
spanning a 05:33:33.7 to 05:37:08.4 and 6 —04:57:57 to —05:51:14. This region is displayed in Fig¬ 
ure During each of the three observing seasons, data were taken about once per night; each 
night’s observations were identical in pointing and filter choice. Because four pointings from the 
four-detector array are required to cover the target field, and we observed at three wavelengths, one 
night’s data comprises 4 x 4 x 3 = 48 image “frames”. As discussed in section 2.2, we have removed 
a number of these frames during processing due to poor data quality. The weather during these 
winters was sometimes poor, so much of our data cleaning strategy revolves around mitigating the 
effect of this bad weather on our light curves by removing suspicious data that are likely to have 
large errors. 


The observations came from UKIRT WFCAM programmes U/06B/H54, U/07B/H60B, and 
U/09A/H33 and were combined into a single WSA combined programme WSERV5. The data 
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in this paper come from the data release WSERV5v20120208, and can be accessed at http;// 
surveys.roe.ac.uk/wsa. 


2.2. Data: Processing and Cleaning 


Before using the images obtained from UKIRT to study variability in the ONC, we needed to 
extract the photometry and then “clean” it to ensure its integrity. The data have been reduced 
and processed for analysis in several stages: first, the raw images in minimum schedulable blocks 
(MSB) were reduced to catalogs and calibrated at the Cambridge Astronomy Survey Unit (CASU; 


Irwin et al. 2004; Hodgkin et al. 2009) Second, the reduced, calibrated images and extracted 


catalogs were transferred to the Wide Field Astronomy Unit in Edinburgh where they were ingested 
into the WFCAM Science Archive (WSA) and multiple MSBs’ data were processed together to 
produce deep images, multi-band catalogs, links from unique sources to epochs, variability tables 
and external surveys (see Hambly et al.||2008 Cross et al.||2009 ) Third, the catalogs were analyzed 
and cleaned after retrieval from WSA. We briefly describe these reduction and processing steps in 
this subsection. 

The data reduction, calibration and archiving that form the first and second step were done 
as part of the VISTA Data Flow System designed to process near-infrared data from WFCAM and 


Emerson et al. 

2004 

), in a similar fashion to the data used in Rice et al. ( 

2012 


Stellar magnitudes are calibrated against local Two Micron All-Sky Survey (2MASS) standard 


sources with high signal-to-noise (see Hodgkin et al. 2009). The WSA produces a catalog J, H, 


K magnitudes and colors for each epoch (see Cross et al.||2009 §2.2 for details), with a maximum 
epoch duration of 720 minutes. Additionally a single master catalog is created from merging J, 
H, K catalogs extracted from deep stacks of all epochs. As a result of the size of this survey, a 
large number of stars have data quality issues (such as saturation, blending, bad pixels, and other 
effects). These stars are given error flags by the WSA pipeline whenever they arise, and we have 
used these error flags to help determine which light curves have a quality high enough to allow 
automatic light curve classification (see Appendix [A]) . 

Some regions in individual nightly catalogs have been removed by WFAU due to obvious 
quality issues. We rejected 137 detector frames (there are 4 in a normal multiframe) because they 
had a dead detector or all the channels were bad, and one channel because it had a value of 0 for 
the seeing, indicating that a non-sensical catalogue was produced. We also rejected 738 detector 
frames when we recalibrated the data (small adjustments to improve the light curves) because the 
adjustments were > 0.05 mag. This is usually an indicator that there is a problem with the frame. 

After obtaining the night-by-night catalog of photometry for all sources, we have applied a 
systematic correction to the pipeline-estimated photometric errors. The correction, derived by 


Hodgkin et al. (2009) and following Rice et al. (2012) appears as follows: 
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where M is the measured (i.e., realistic) total error, E is the pipeline-estimated photometric error, 
and c = 1.082 and s = 0.021; c and s are empirically determined. The observational error as a 
function of magnitude is plotted in Figure for each band. 


We removed individual exposures showing unusually high errors in photometry relative to 
typical observations^ To quantify the fidelity of a given exposure, we have selected a subset of 
likely constant “reference” stars using a version of the Stetson index (described in 13.1 as part of 
our variability criteria); reference stars were also required to be between magnitudes 13 and 15.5, 
in order to guarantee high signal-to-noise while widely avoiding saturation effects. An example 
reference star light curve is shown in Figure At least 20 reference stars were selected in each tile 
to ensure that no single star could disqualify an exposure. These stars’ median observed magnitude 
was tabulated, and then each night’s photometry was compared to this median magnitude for each 
night. We have assessed photometric fidelity on each night by measuring the fraction of reference 
stars that fall within 0.05 mag of their median magnitude. In the vast majority of observations, 
fewer than 10% of reference stars (i.e., > 90% fidelity) deviate by 0.05 mag; we have chosen to 
discard all exposures in which more than 20% of reference stars deviate from their typical magnitude 
by at least 0.05 mag (i.e., fidelity less than 80%) as unreliable for our analysis. 


Next, after noting that most “clean” light curves still had one or a few flagged datapoints in 
each band, we removed those minority errant datapoints (in this context, data “flags” note data 
quality issues such as deblending, edge effects, saturation, diffraction spikes, etc.). This procedure 
analyzed each band (of J, H, and K) independently: for each band on each star, the number of 
flagged versus flag-free datapoints is counted. If a star has > 90% flag-free data in a given band, 
then any flagged datapoints are removed. Finally, all datapoints with error > 0.5 mag are removed, 
because such data occasionally clutter light curves that otherwise contain good data. 


From the above procedure, we extract 14,728 lightcurves meeting the minimum quality re¬ 
quirement {Q = 0) or better; a distribution of their median magnitudes is shown in Figure for 
each band. Example light curves of various types can be seen in Figure 


3. VARIABLE SELECTION 


In this study we focus on stars with significant variability. We have chosen to use the Stetson 
index (Stetson 1996 and citations therein) to quantify variability seen in our JHK light curves; 
this index is ideal for measuring variability that is correlated between bands, and so it makes the 


^ In this context, an “exposure” refers to photometry at one band that was collected using one pointing on one 
night; a typical exposure covers four tiles/“frames”, spaced as in Figure]^ though some exposures had one or more 
frames removed by WFAU as described above. Thus, a typical night had 4 x 3 = 12 exposures, for 4 pointings and 3 
bands. 
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best use of our three-band datsj^ In order to characterize the periodicity of these stars’ variability, 
we have used the Lomb-Scargle Periodogram (see Press Sz Rybicki 1989 for our implementation), 
which searches for sinusoidal periods in irregularly time-sampled data, as well as the Fast Chi- 
Squared (fy^) method developed by Palmer (2009) in specific cases where it was clear that the 
variability was not sinusoidal (e.g., eclipsing binaries). 


3.1. Stetson Index 

The Stetson variability inde:x ^[Stetson (1996) is defined as follows: 


5 = 


Ef=i 9 i sgn (Pj) ^/\W\ 

E n 5 

i=l 9i 


( 2 ) 


where p is the number of pairs of simultaneous observations of a star, and Pi = is the 

product of normalized residuals of two observations (where the normalized residual is defined 


5i = 



( 3 ) 


for a given band; n is the number of observations in that band, m, represents the magnitude at a 
given epoch, ai is the observational uncertainty on rrii, and m is the average magnitude of the star). 
Most stars in our analysis have many simultaneous J, H and K measurements. In the Stetson index, 
three simultaneous observations are treated with = | for each pair of measurements, giving a 
total weight 2 for the set of three pairs. Many stars lack complete data in three bands (for example, 
they are too faint at J or too bright at K)] for these stars, a two-band or one-band Stetson index 
has been calculated where appropriate, with weights set at gi = 1. 


The Stetson index is designed to identify significant variables at S' > 1.0 when noise properties 
of data are well characterized; truly nonvariable stars will cluster around 5 = 0. In practice, 
previous authors have found selection criteria of 5 > 0.5 (Carpenter et al.||2001 Morales-Calderon 


et al. 2011; Flaherty et al. 2013) to be useful for determining minimal variability criterigj^ 


Because of the large range of photometric quality issues that arise in a major dataset such as 
ours, we have categorized all our sources into one of three quality “bins”. Briefly, the stars are 


^ While the Stetson index de-emphasizes anti-correlated variability in principle, we find that stars showing obvions 
anticorrelated events have generally correlated variability in practice, leading to large Stetson indices; thus, the use 
of Stetson is useful even for stars with some anticorrelated variability. See Appendixfor further discussion. 

^We rename the J Stetson index to “S" to avoid confusion with the J wavelength band, following previous authors’ 
usage. 

'^Low-amplitude variable candidates are discussed in Appendix B 



























divided into Q = 2, Q = 1, or Q = 0. Q = 2 stars have completely clean data in all three bands; 
Q = I stars have one or two bands with completely clean data; Q = 0 stars have some quality issue 
in all three bands. 

Each quality bin has different variability selection criteria, due to the varying confidence in 
light curve fidelity. In all cases, the minimum variable criteria are S' > 1.0 and at least 50 detections 
in at least one band; stars in lower “quality” bins were required to meet further criteria. Q = 2 
and Q = 1 stars have been extracted algorithmically, while Q = 0 stars required visual light curve 
inspection to be included in our list of variables. Full details are provided in Appendix [A} We were 
very conservative in selecting Q = 0 variables, taking care to only include light curves that are 
clearly not false signatures from instrumental or weather-related effects. No Q = 0 stars (variable 
or non-variable) are included in any statistical analyses presented later in this paper. 

Using the Stetson Index and the selection criteria described in Appendix [Aj we identify 1203 
significant variable stars, out of 14,728 stars tested. The breakdown of these 1203 variables into 
quality classes is described in Table as Q = 2, Q = 1 and Q = 0. The distribution of Stetson 
indices among the Q = 2 sources is clustered around a median index S = 0.487 as shown in Figure 

El 


3.2. Periodicity Analysis: LSP, 


Much observed variability among stars is periodic. To detect periodic variability in our data, 
we have applied the Lomb-Scargle periodogram to the J, H, and K band light curve of every single 
variable candidate. The Lomb-Scargle periodogram was developed to identify and characterize 
periodic, sinusoidal signals in irregularly time-sampled data, and has been used widely in other 
studies of variability in young stars (see, e.g.. 


Attridge & Herbst||1992 

Stassun et al. 

1999 

Rebull 


2001 ). 


Among light curves that suggest a periodic (often high-frequency) non-sinusoidal signal, we 
have used the Fast Chi-Squared (fy^) method developed by Palmer (2009) to get a more accurate 
period fit than the Lomb-Scargle periodogram can provide. This is most effective for eclipsing 
binaries, W UMa “contact binary” variables, as well as stars that show any sort of abrupt “dip” or 
“burst” that is periodic but surrounded by a flat, “quiescent” light curve. 


3.2.1. Automatic period detection algorithm 

In order to efficiently detect strong periodic signals in our large dataset of light curves, we have 
created a filter to automatically classify stars as periodic or non-periodic based on their Lomb- 
Scargle periodogram strength. The Lomb-Scargle periodogram returns a plot of strength versus 
frequency; periodogram strength can be transformed into a corresponding “false-alarm probability,” 
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assuming Gaussian errors around perfect sinusoidal signals. Example periodograms for two stars 
are shown in Figure 

We have calculated Lomb-Scargle periodograms in the range 2 — 50 days for all candidate 
variable stars; if a star showed periodogram strength exceeding 15 in a single band (i.e., false 
alarm probability, or FAP, < 0.0002) or exceeding 12 at the same period, within a tolerance of 
5%, in two or more bands (i.e., FAP < 0.0041 in each band independently), then we classified it 
as either a “periodic variable” (if it was in Q = 2 or Q = 1 and was selected as a variable) or 
a “subjective periodic candidate” (if it was in Q = 0 and was selected as a candidate variable, 
pre-inspection). These filters have been tuned by hand after much trial-and-error and inspection of 
individual light curves, and are meant to detect the strongest, repeating periodic variables rather 
than being sensitive to every statistically-significant periodic signal. Because of the large sample 
size of our dataset, this optimization is preferable to imposing a target false-alarm probability from 
the top down. 

The lower limit of 2 day periods comes from a Nyquist limit (i.e., a fundamental sampling 
limit) on our ability to identify short-period sinusoidal signals in our data; excluding period fits 
below 2 days guarantees that true periods above 2 days are not aliased to below 2 days, at the 
potential cost of aliasing some (poorly-sampled) sub-2-day periods upwards. This aliasing effect 
arises from the fact that a given set of discretely sampled points can typically be fit by two different 
sine waves: one with a frequency above the Nyquist frequency of half the sampling frequency, and 
one with a frequency below the Nyquist frequency. The upper period limit of 50 days was chosen to 
ensure that periodic variability was reliably repeatable, and would be observed to repeat multiple 
times in the data. 


3.3. Results of Variability Selection 


The breakdown of the 1203 variables into quality bins and periodic/nonperiodic bins is shown 
in Table In this section, we briefly characterize the overall variability and periodicity properties 
of our variables, as well as compare our set of variable stars to a number of catalogs from the 
literature, allowing us to comment on the accuracy of our period search, identify cross-matched 
identifications in SIMBAD and other datasets, and assign Spitzer photometry and evolutionary 
classes to 90% of our stars. The coordinates, SIMBAD cross-matched names, median WFCAM 
photometry, Spitzer photometry, and evolutionary classes of all 1203 variables are listed in Table 
Stellar classes are taken from Megeath et al. (2012) and are used extensively in our analysis in Q 
SED classification suffers from a degeneracy between the extinction and source inclination (both 
poorly constrained), which in particular affects the determination of source luminosity ( ]Fischer 
et al.|[20l2) ) and possibly also the YSO spectral class. This should be kept in mind in the following 
discussion. 
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History of infrared monitoring campaigns in the ONC 


Carpenter et al. 2001 

T This paper 
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Fig. 1.— A history of infrared monitoring campaigns in the Orion Nebula Cluster published to 
date, with each band (J, H, K, [3.6], [4.5]) shown in a different color. The survey presented in this 
paper lasted from October 2006 to April 2009. 


Notes- “ YSOVAR refers to Morales-Calderon et al. (2011). 


Table 1. Breakdown of variables into period & quality bins 



Periodic 

Non-Periodic 

P-hNP 

Non-variable 

Total 

Q=2 

231 

318 

549 

1799 

2348 

Q=1 

116 

192 

308 

936 

1244 

Q=0 

193 

153 

346 

10790“ 

11136 

total 

540 

663 

1203 

13525 

14728 


Note. — “ For Q = 0, the “Non-variable” column enumerates 
stars not selected as variable, but which are not necessarily non¬ 
variable; see Appendix [A[ 
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Fig. 2.— An illustration of the observed tiles in our field, with the effective number of observations 
in each band written in each tile, i.e., the number of observations that remain after the data cleaning 
process described in §2.2[ The four gray tiles illustrate a single detector “pawprint” of WFCAM; to 
image the entire square degree, four pointings of this pawprint are required. The spacings between 
WFCAM chips cause some overlap in the pawprints, so a small area of the field is imaged multiple 
times per night. 
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Fig. 3.— Spatial distribution of the 1203 variable stars selected in ^ Each red marker on this map 
shows the position of one variable. Our observations are centered on the Orion Nebula Cluster and 
cover a ~ 0.9 x 0.9 deg^ field. The brightest stars in the ONC, including the well-known Trapezium 
stars at the center, are saturated in our images, preventing us from studying their variability. The 
highest density of variables runs along the north-south ridge corresponding to a filament of active 
star formation. For comparison, a detailed map of dust optical depth in this region can be found 
Lombardi et al. (2014). Credit for underlying M42 near-infrared {Z, J, Kg) image: ESO/J. 
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Fig. 4.— An example light curve of a constant (reference) star, to demonstrate the level of 
variability detectable in this study. Errors in each band are about 0.02 mag, depicted by whiskers 
around the photometry points. The vertical dashed lines mark where “empty” stretches of time 
were omitted from the plot. In the bottom-right panel, a color-color diagram is shown, with the 
solid lines denoting a main sequence locus (Koornneef||1983 thin) and a CTTS locus (Meyer et al. 


1997, thick), and dashed lines showing reddening vectors using the ISM reddening law of Rieke & 


Lebofsky (1985), extending from the extrema of the main sequence and CTTS curves. These solid, 


dotted, and dashed lines appear in all light curves and color-color diagrams throughout this paper. 
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Fig. 5.— Eight example light curves highlighting the variety of periodic and nonperiodic behaviors 
seen in our data. The top four panels show periodic variables folded on their own periods; the 
bottom four panels show nonperiodic variables displayed without folding. 
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Fig. 6.— Observed rms error as a function of magnitude, measured for nonvariable stars. The 
photometry has roughly 2% error down to a limit of J = 16, H = 15.3, and K = 14.8, and better 
than 5% error to J = 17.3, H = 16.7, and K = 16.2. We restricted our highest-quality subset of 
variable stars (“Q = 2”) to data within 11 < J < 17, 11 < 77 < 16 and 11 < 77 < 16, among other 
criteria. 
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Fig. 7.— Distribution of mean magnitudes among sources, per band. The empty (light gray) 
region shows the distribution for all sources detected at least 50 times; the filled regions show only 
stars that meet the stricter Q = 1 ot Q = 2 quality requirements. 
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Fig. 8.— Distribution of the three-band Stetson index versus median H magnitude, for all 2,348 
Q = 2 stars. The distribution peaks at a median value of S' = 0.487; we select S > 1.0 as our 
strong variable star definition. Stars between 1.0 < S < 0.55 are excluded from our analysis but 
are discussed briefly in Appendix 150 stars in this sample have S > 5 and are not shown in this 
plot. 
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Fig. 9.— Three-band Lomb-Scargle periodograms for two stars, plotted in period space from 2 to 
100 days. From top to bottom: J (blue), H (green), K (red). Left column: periodogram of star 
ONCvar 40, a rotating star (light curve shown in Fig. 14). It shows a strong, sharp peak at 16.6 
days, indicating a clean periodic signal. Right column: periodogram of star ONCvar 553 (light 
curve shown in Fig. 30). Its 11.0 day period is significant, but has notable structure, indicating 
that the lightcurve has substantial irregular variability superposed on a periodically varying signal. 





























Table 2. Basic Properties of 1203 Significantly Variable Stars 
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Note. — Sample; the remainder of the table is available electronically. Column Q refers to quality class, and column P denotes whether a star is 
periodic (1) or not (0). J, H, K values listed are median values, with representative errors on the median measurement. ^ P=“Protostar”, D=“Disked”, 
ND= “Non-disked” (classes derived from|Megeath et al.|l|2Ql^). 
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3.3.1. Global Variability Properties 


By considering the total number of stars (variable and nonvariable) that met the Q = 2 or 
Q = 1 criteria, we can study the fraction of stars in our input catalog that are variable, periodic 
or not. This approach minimizes selection biases, because all Q = 2 and <5 = 1 stars have been 
subject to nearly the same data quality, magnitude, PSF, and spatial cuts, the only difference being 
the number of bands which passed those tests. Results are statistically indistinguishable when the 
Q = 1 and Q = 2 groups were considered separately, so for global considerations we combine them 
while excluding Q = 0 data. There are 3592 Q = 1&:2 stars of which 857 are variables with 5 > 1. 
Overall, the fraction of our 3592 stars that are variable is thus 24.2%. Of these, the fraction of 
variables that we found periodic is 0.41. Therefore, the percentage of stars in the overall dataset 
that are periodic variables is 9.9%. 


Our results can be compared to other NIR monitoring campaigns of this region: Carpenter 


et al. (2001), who covered a much larger area (0?8x6°, oriented north-south) at JHKg to a shallower 


magnitude of Kg = 14.8; and Morales-Calderon et al. (2011), who studied a similarly-sized field at 
[3.6] and [4.5] Spitzer bands that overlaps with most of our observed field. Carpenter et al. (2001) 
found 233 periodic variables among 1235 variables in their 17,808 star sample (i.e., ~ 7% of their 
overall sample was variable, a fraction ~ 0.19 of their variables were periodic, giving ~ 1.3% of 
the overall star sample as periodic); [Morales-Calderon et ah (2011) identihed a higher proportion, 
501 periodic variables among 1259 variables in a 2069-star sample (i.e., ~ 60% of the pre-selected 
sample were variable, a fraction ~ 0.40 of their variables were periodic, giving ~ 24.2% of the 
pre-selected star sample as periodic). The clear differences between the three surveys are due to 
several factors: a shallower survey will probe only the upper, less-populated part of the initial mass 
function and will not penetrate as deep into higher extinction regions, while a survey at longer 
wavelengths is likely to target younger sources that are more likely to be variable. Beyond this, the 
Morales-Calderon survey focused on a pre-selected star sample. 


Because our variability and periodicity cuts are conservative and specifically exclude very low 
amplitude variability and short periods, these fractions can be taken as lower limits for the intrinsic 
periodicity and variability fractions of stars in our field. 


3.3.2. Literature Cross-Matching 


Because the ONC is an extremely well-studied region, we have been able to access numerous 
catalogs containing previously-published data relating to the variable stars we present in this paper. 
We have obtained previously measured periods, Spitzer-deVwed evolutionary classes, and XMM- 
Newton-derived evolutionary classes, from the literature. 


We compare the periods found via our technique described in [ 3.2.1 to six catalogs of periods 
from the literature: Carpenter et al. (2001), Morales-Calderon et al. (2011), Herbst et al. (2002), 































- 21 - 


Parihar et al. (2009), Rodn'guez-Ledesma et al. (2009), and the General Catalog of Variable Stars 


( ]Samus et al.|[2009 ). The Herbst, Parihar, and Rodriguez-Ledesma campaigns were all done in the 
optical, finding periods spanning | to 22.2 days, 0.9 to 19.1 days, and 0.09 to 16.9 days, respectively. 
The Carpenter study was a NIR (J, H, Kg) project, while Morales-Calderon used the [3.6] and 
[4.5] Spitzer IRAC bands. 

Of the 539 periodic variables we discuss in this paper, a comparison with the literature shows 
that 297 of these stars are already known to be periodic, while 242 are previously unknown periodic 
stars. Where our period data overlap with the literature, they largely agree, as shown in Figure pR) 
There is disagreement where the literature periods are shorter than two days because our sampling 
is not designed to be sensitive to these shorter periods. For overlapping stars, the fraction of our 
periods that agree with published periods within ±10% is tabulated in Table The agreement 


is greatest for the two infrared studies (Carpenter et al. 2001 and Morales-Calderon et al. 2011); 


this is indicative that infrared and optical variability do not systematically arise from the same 
mechanism(s). 

The overlap is not perfect: 194 of our non-periodic variables were previously reported in the 
literature as periodic. Because the bulk of previous period studies in the ONC were done at optical 
wavelengths, and some of them had far finer time resolution than our ~nightly survey, the non¬ 
recovery of these 194 periods is likely attributable to aj the observation that infrared variability 
is often more irregular (or simply different) than optical variability (see simultaneous optical and 
infrared monitoring of YSOs in Cody et al. |2014 ), due to the emission of the disk at infrared 
wavelengths, and b) our survey is not designed to recover periods below two days. Indeed, among 
these 194 stars, 44 have literature periods below two days (22% of the sample), so case b) accounts 
for almost one-fourth of the missing periodic stars. 


In order to determine evolutionary classes for our variable stars, we made use of the Megeath] 


et al. (2012) study. This study identifies YSOs based on their infrared data from 2MASS (J, 


Table 3. Period agreement for each catalog 


Period catalog 

Fraction that agree 
within ±10% 

Number of matches 
to this catalog 

GCVS 

89.58% 

of 48 

Carpenter±01 

94.59% 

of 74 

Morales-Calder6n± 11 

90.00% 

of 40 

Herbst ±02 

79.74% 

of 153 

Parihar±09 

77.05% 

of 61 

Rodriguez-Ledesma±09 

66.67% 

of 96 
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H, Kg near-infrared bands) and Spitzer (IRAC [3.6], [4.5], [5.8], [8.0] and MIPS [24] mid-infrared 
bands). Stars with significant infrared excesses are identified as YSOs, and those meeting a spectral 
index criterion of a > —0.3, where a = dXFx/dX, are further identihed as protostellar YSOs. 
The subsequent catalog classihes all sources as “Protostar”, “Disked”, “non-Disked”, or “data too 
imprecise to determine”. As noted in the first paragraph of Section 3.3, this class determination 


suffers from inclination degeneracy. We have cross-matched our variable stars to this Megeath et al 


(2012) catalog and found that 90% have a corresponding evolutionary class. Of our 1203 variables. 


47 match to protostars, 641 to disked stars, and 390 to diskless stars, with the remaining 125 either 
matching to stars with Spitzer colors too imprecise to determine a class, or failing to match to a 
Spitzer source altogethei|^ A comparison of periodicity (among Q = 2 variables) to evolutionary 
class shows 27% of protostars, 34% of disked stars, and 55% of diskless stars as periodic in our 
data. 


We have also cross-matched our sources to the X-ray identified YSOs in Pillitteri et al. (2013), 


which was an XMM-Newton study of the ONC. This catalog includes X-ray identified YSOs, and 
assigned them classes (“Protostar”, “disked”, and “Class IIP”) based on the Megeath et al. (2012) 


classification; through this, young stars can be identified independent of the presence of a disk. 
Class III status was assigned to X-ray sources with no IR-detected disk. Of the 850 sources in this 
X-ray catalog, 270 matched to stars in our study, of which 166 had estimated evolutionary classes. 


4. Analysis and Results 


In addition to the overall variability statistics presented in the previous section, we have an¬ 
alyzed these stars’ variability in light of prior information about the stars themselves as well as 


near-infrared variability mechanisms. As discussed in [3.3.2 the Megeath et al. (2012) catalog 


provides Spitzer-derived evolutionary classes for most of our variables; we analyze the dependence 
of variability on evolutionary class in §4.H and the dependence of period distribution on class in 


[4.2 Carpenter et al. (2001) found that a star’s slope in the H — K,K color-magnitude diagram can 
be used to constrain or identify the physical mechanism causing its variability; in §4.3] we analyze 
the color slope behavior of the bulk population as well as individual stars. 

To give some context about the physical properties of such young stars and their disks, a 


reference such as Eisner et al. (2009) can be consulted. The PMS stars RY Tau, DC Tau, and RW 


Aur A were measured to have inner disk walls in the range 0.14 — 0.18 AU for these three systems, 
with disk wall temperatures of 1260 — 1850 K. Such systems have orbital periods of 5 — 20 days for 
their inner disk walls; other typical, lower-mass stars have inner walls at 0.1 AU, with a roughly 
11 day orbital period. 


®These non-matches did not have a counterpart in the Megeath catalog within 1 arcsecond. Most of the non¬ 
matches are in crowded, easily confused regions, and the UKIRT seeing is typically better than Spitzer’s images, so 
we are not concerned that we have included spurious sources in our variability catalog. 
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Fig. 10.— Comparison of previously-published periods with periods derived in this study, for the 
297 stars in which the data overlap. A strong general agreement is seen between literature periods 
and our periods, except where the literature periods are shorter than two days; such periods get 
aliased to longer periods in our data. The two nearly horizontal curves that appear below two 
days on the y-axis are due to such aliasing, which is expected given the time resolution of our 
study; a black dashed line drawn at y = 2 is drawn to indicate that periods below this are affected 
by aliasing. In addition to the overlapping periods, we hnd 242 periodic variables unknown to 
have periods in the literature. We expect that the accuracy of these newly reported periods is 
comparable to that seen among our periods for the previously-known periodic stars. 








- 24 


A color-color diagram showing the median J — H^H — K values of our variables, by evolu¬ 
tionary class, shows that disked stars tend to have larger near-infrared excesses than diskless stars 
(Figure [IT|). Note that such a color-color diagram only captures the median color of a star in our 
observations; in reality, individual stars move substantially throughout color-color space, making 
near-infrared excesses variable (Rice et al.||2012). 


4.1. Variability as a Function of Class 


Young stars are known to be more variable than main-sequence stars (see e.g., Morales- 


Calderon et al. 2011; Herbst et al. 2002); using the evolutionary classes associated with the stars 


in our sample, we can refine this observed trend and investigate the variability associated with 
each class. As discussed in 13.3.2 we obtain evolutionary classes for 90% of our variables from the 


Megeath et al. (2012) catalog; many non-variables in our sample also have classes. A large fraction 


of these are identified as protostars or disked stars: among 2373 Q = 1 + 2 stars (both variable 
and nonvariable) with Megeath identihcation, 27 are protostars and 639 are disked, while 1707 are 
nondisked. The Stetson index statistics presented below are from this sample. 


The disked sources have a median Stetson index of 3.0 =fc 2.3 p| (Figure [I^ , while the protostars 
are fewer in number and have a median value of 5 = 7.6 ± 4.8. Since the Stetson index is an 
aggregate of behavior across all bands, we analyze variability strength at different wavelengths 
among variable stars by inspecting amplitudes in each bancQ In the following analysis we restrict 
ourselves to the 522 Q = 2 variables with Megeath classes, of which 11 are protostars, 288 are 
disked, and 223 are nondisked. Again we find the protostars are far more variable: median AK 
{Kmax — Kmin) for the Variable protostars is 0.61 ± 0.34 while the corresponding value for the 
other disked variables is 0.21± 0.10. Non-disked stars are typically much less variable: median 
AK = 0.09 ± 0.02. If the Q = 1 protostars are included alongside the Q = 2 protostars, giving 27 
protostars in this analysis, the median AK {K^ax — Kmin is 0.492 ± 0.29, which is consistent with 
the Q = 2-only case. 


The color terms A(J — H) and A{H — K), measured as the difference between maximum and 
minimum color values after iteratively clipping 3a outliers, show similar trends: color shifts are 
smaller for more evolved stars. As an example we show the results for changes of the H — K color 
in Figure 13 We primarily use the A{H — K) color term to probe changes in the disk excess 
Indeed, we find the median A{H — K) = 0.26 ± 0.13 for the protostars, A{H — K) = 0.12 ± 0.04 


®The presented uncertainties on all median values are calculated as the median absolute deviation. 

^ When calculating Amag for a given star, extreme (3 (t) outliers in a star’s lightcurve are iteratively clipped so 
that errant data cannot bias the magnitude ranges upwards for low-amplitude variables. 

^ J — K is equally as useful as H — K, but as the redder stars in our sample are very faint and noisy at J, the 
H — K index is more reliable. 
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All variables, colored by Megeath+2012 class 



median H-K 


Fig. 11.— Median J — H,H — K color-color diagram of the 1,077 variable stars discussed in 
this paper that have an evolutionary class from Megeath et al. (2012). As in Figure [^s color-color 
diagram, the solid lines in the bottom-left show the main sequence ( Koornneefjl983| ) and the CTTS 
locus of Meyer et al. (1997), while dashed lines extending to the upper right show reddening vectors 
of typical ISM dust (Rieke & Lebofsky 1985). Stars with disks tend to have larger near-infrared 
excesses than diskless stars, with protostars showing the largest excesses. Note that only the median 
values of J — 77 and H — K are plotted for each star; in reality many stars change their position in 
color-color space, with less-evolved stars showing the greatest color changes (see Figure [l^. 
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0 1 2 3 4 5 



0 1 2 3 4 5 


Stetson Index 


Fig. 12.— Stetson index for each class. Median Stetson indices are computed for the aggregate 
Q=1 + Q=2. The x-axis is limited to the range 0 — 5 to highlight the differences between the 
distributions; in the top panel, 16 sources lie beyond 5 = 5; in the middle panel, 232 sources he 
beyond 5 = 5; in the bottom panel, 22 sources lie beyond 5 = 5. 




- 27 - 


for the disked variables and A{H — K) = 0.06 ±0.01 for the diskless variables. 


The Spitzer sample is able to identify diskless objects, but it is limited by its inability to dis¬ 
tinguish diskless YSOs from stars that have evolved onto the main sequence. To better understand 
the variability properties of young diskless stars without contamination by main-sequence stars, 
we use an XMM-Newton X-ray selection, as X-ray data provide a sample of YSOs unbiased by 
the presence of a disk. The central portion of the monitored field was subjected to three XMM 
pointings averaging about 50 ks each. X-ray detections of star-like sources with IRAC [3.6]-[4.5] 
color < 0.3 are likely to be Class III YSOs ( Pillitteri et al.|[2013 ® so they are identified as such in 


their XMM-Newton-selected catalog. 270 XMM X-ray sources were identified with Q = 1 ± 2 stars 
in our monitored held; 166 of these sources are identihed as young stars in the XMM catalog. This 
includes three Class I objects, and a roughly even split between the number of Class II {N = 90) 
and Class III {N = 73) objects detected in X-rays. Statistics below are limited to the Q = 2 stars 
in these samples: N=28 (Class II) and N=49 (Class III). The median Stetson indices and AK 
values observed for the Class II stars in this sub-sample are consistent with the Spitzer-selected 
disk sample {S = 2.54 ± 1.74 and AK = 0.20 ±0.10) conhrming that the X-ray sample is relatively 
unbiased with respect to the NIR variability. Meanwhile, the disk-free young stellar objects (“Class 
III” ) are far less variable by both criteria (median S = 1.14 ± 0.62 and median AK = 0.09 ± 0.03). 
The median AK = 0.09 value is fully consistent with variability arising only from starspots (e.g.. 
Carpenter et al.|2001 Wolk et al.|2013 ). In both the X-ray sample and the Spitzer-selected sample, 
we hnd that less-evolved objects are by far the most variable. 

It is further interesting that protostars are significantly more variable than evolved disked sys¬ 
tems. The distinguishing feature that separates protostars from evolved disked stars is the presence 
of a dense outer circumstellar envelope, which near-infrared observations at K band (2.2 /rm) are 
not sensitive to. The JHK bands should only be measuring the stellar photosphere, the star-disk 
interface, and the innermost edge of the dust disk. Thus, the inner disk edge and/or the star-disk 
interface must have a fundamentally different structure between protostars and disked stars, with 
the protostar inner environment more unstable, leading to either accretion/disk geometry events, 
dust-occulting events, or (most likely) a combination of both. The instability is strongest for the 
least evolved systems. Note that, by definition, protostars are in the evolutionary phase where 
they accrete most of their mass; it is therefore expected that the least evolved young stellar objects 
show stronger and likely more frequent instabilities and variability due to the accretion process. 
An alternative interpretation of this result is that the distribution and stability of the hottest cir¬ 
cumstellar dust graim , thought to be responsible for A'-band excesses, is different for protostars 


when compared to typical disked objects. 

In Tables and we show the variability statistics, including derived color slopes, for all 


®We do not use the deeper Chandra data available because the Chandra-observed region is much smaller than the 
XMM-observed region and hence more subject to spatial bias. 

^°The solids at the inner edge are refractory oxide grains (e.g., CaO, AI2O3, Ti203, MgO) 


(Pollack et al. 


1994 
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periodic and nonperiodic stars in our dataset, respectively. We sort each table by Megeath et al. 
(2012) class. 


4.2. Period distributions by Class 


Periodic behavior is frequently found in variability studies of YSOs, and searches for periodicity 
have been used, e.g., to study the evolution of angular momentum loss throughout early stellar 
evolution. Periods are regularly found for both disked and diskless objects; well-sampled surveys of 
some star-forming regions have found a bimodal distribution of periods in which disked stars have 
periods of ~ 8 days, while mostly disk-less stars have periods under 4 days and diskless stars tend 
to have shorter periods overall^ In the near-IR, the first signihcant study of periodic behavior 
was that of Carpenter et al.| ( |2001[ ). They found a more or less continuous distribution of stellar 
periods spread between 2 and 10 days. The distribution peaks with about 45 stars having periods 
between 2-3 days, and falls smoothly to 15 stars with periods between 10-11 days. The peak value 
is probably augmented by stars with true periods shorter than the sampling frequency, which would 
imply a somewhat flatter distribution. The bulk of the periodic stars in that sample did not have 
K band excesses indicative of disks. 

In Figure seven examples of cleanly rotating spotted stars are shown. 

In Figure [Tsl we show the period distribution for 471 periods found using the Lomb-Scargle 


method as described in 13.2.1 All periodic sources (regardless of quality status) with a Spitzer- 
derived evolutionary class are included in this histogram, provided they have periods below 20 days. 
While the distribution for the protostars is poorly sampled, the disked and diskless distributions 
are well populated and intriguing in a couple of respects. First of all, the disked and disk-free 
populations are distinct; the disk-free population has more short periods (below 5 days) with 
respect to longer periods than the disked population. The distribution seen among the non-disked 


sample is very similar to that found by Carpenter et al. (2001). Here, we find about 5 times more 


stars with periods in the 2-3 day range as in the 10-15 day range. One caveat about the period 


distribution we hnd, which is common to Carpenter et al. (2001), is that some periods are likely to 


have been moved from below 2 days to above 2 days via aliasing. 

The distribution of periods for disked stars is significantly different from that of the diskless 
stars: a two-sided Kolmogorov-Smirnov test gives a near-zero probability that these two groups of 
periods were drawn from the same distribution. A small enhancement of disked objects at periods 
between 5 — 10 days marks a key difference between the disked and diskless distributions; there is 


^This bimodality has been hailed as evidence of magnetic coupling between the star and its disk (e.g., 


Attridge & 


Herbst|1992 Edwards et al. 1992, Rebullj2001 ). In this model, the star-disk coupling prevents the star from spinning 


up while the radius shrinks during the contraction phase of the Hayashi track ( Hayashi &: Nakano||1965 1. Once the 
disk dissipates, the stellar magnetic field is decoupled and the star is free to spin up as mass settles to the center of 
the star. 





















Table 4. Variability properties of periodic stars 
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stars fall on overlap regions and are observed twice per night. m(x, y) denotes the slope in the specified color-color or color-magnitude diagram; a 
slope is only included if the orthogonal distance regression fit yielded a confidence of better than 10% on the fit. S denotes the Stetson Variability 
Index. 












Table 5. Variability properties of nonperiodic stars 
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K mag range for Q=2 variables 


2.0 



0.5 1.0 1.5 2.0 



0 ^ ---^^- 

0.0 0.5 1.0 1.5 2.0 


AK mag (outlier-proof) 


A{H-K) color (outlier-proof) 


Fig. 13.— Distributions of variability amplitude at K band and H — K color for protostars, disked 
stars, and diskless stars. Typical amplitudes are much higher for less-evolved stars versus more- 
evolved stars. The median K amplitudes for each class are 0.61, 0.21, and 0.09, for protostars, 
disked stars, and diskless stars, respectively; median H — K amplitudes are 0.26, 0.12, and 0.06 
respectively. 
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Period-folded light curve 


bght curve (MJO-54034) 




J-H VS. H-K 

1.0 

00 

ot 

0 7 

04 

✓ 

✓ 

✓ 

y 

y 

<- 


0200250 MO.HO 40 

l.M 

l.OO 

i.es 

l.M 

i.n 

1.70 

1.6S 

✓ 

y 

■ ^ 

. / ■ 


0.4S 0.00 005 1 00 

lO 

00 


04 

0 7 

04 

✓ 

✓ ** 

^ < 


0200250 MD.MO 40 

lO 

00 

04 

07 

04 

y 

y 


0200250 MO.HO 40 

lO 

00 

04 

0 7 

04 

y^^ 

y 

y 

y 

c . , y, - 


0200250 MO.HO 40 

lO 

00 

04 

07 

04 

✓ - 
y 
y 
y 

y ^^000^ 

<■ 


0200250 MO.MO 40 

lO 

00 

04 

- ' ' y- 

y 

y 

y 

07 

04 



Oi2oa]ai MO.HO « 


K VS. H-K 



0 iro uaiui.MOJiDM 


Fig. 14.— Light curves of seven stars that show clean, rotationally modulated light curves 
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no peak in the latter half of the distribution as seen in the I band studies (e.g., Herbst et al.||20()2 


Stassun et al.||2000). However, one difference seen between our disked sample and previous studies 


of disked star periods is that fully 1/4 of our disked sample have rotation periods between 2 — 5 
days. In contrast, earlier studies found that disked stars (usually identified via H — K excess) tend 
to be less likely to have periods below 5 days (e.g., Herbst et al. 2000 Rebull||2001 ) 

On the other end of the period distribution, stars with periods longer than 15 days are seen 
in both the disked and non-disked populations. The diskless long-period stars in this sample could 
just be slow rotators. But among disked stars, the longest periods are up to 90 days (see 16.3 e.g.. 
Fig. 36) and are accompanied by color changes consistent with changes in the disk or dust along 


the line of sight. These long-period, color-changing light curves cannot be modulated by the star’s 
rotation period unless the object rotates very slowly - instead, they trace disk phenomena with 
relatively long orbital periods, such as orbital variation in disk extinction or dynamically modulated 
accretion. 


4.3. Color slopes constrain physical mechanisms 


By studying the trajectories traced out by variable stars through color-color and color-magnitude 
space, we are able to add additional constraints to the physical mechanisms that cause their vari¬ 
ability (Figure [Tb]). Accretion disk models generally indicate that motion up and to the right in 
a H — K, K or J, J — H color-magnitude diagram is a result of the additional flux re-radiated at 
long wavelengths as the total flux increases due to accretion, whereas motion down and to the 
right is a result of dust reddening. In the JHK color-color diagram, the dust reddening tracks 


were well established by Lada & Adams (1992) while a locus occupied by classical T Tauri stars 


was identified by Meyer et al. (1997). In the color-color diagram, unlike in the color-magnitude 


diagrams, the CTTS locus and reddening vectors are not orthogonal, and many young stars have 
color-color diagram trajectories with slopes between the CTTS locus and the reddening vector, as 
seen in Cyg OB7 (Rice et al. 2012). Illustrations of the different color-dependent effects on the 


spectral energy distributions (SEDs) of young stars can be found in (Robitaille et al. 2007). 


Carpenter et al. (2001) were the first to analyze NIR color space trajectories in young variable 


stars. They used this approach as the basis for their findings that most stars’ variability is consistent 
with starspots and that stars with high color variability tend to show changes in dust-reddening 
(slopes between 30 and 90 degrees), with a small minority showing changes in disk geometry or 
accretion rate (slopes between -90 and -45 degrees). In the current study, we have examined the 
color trajectories over our longer-baseline data. 

To determine the color “slope”, we use the orthogonal distance regression package included 
in SciPy (scipy.odr), which fits a line to values that have comparable uncertainties in both the 
color and magnitude (or second color). The ODR code returns both a slope and an esti mated error 


on that slope. We calculate such slopes for • Following 
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Fig. 15.— The distribution of periods below 20 days for protostars (top panel, green), disked 
stars (middle panel, red) and non-disked stars (bottom panel, blue). A modest secondary peak in 
the disked periods is seen at ~ 8 days; such a bimodal distribution has been observed before for 
optical periods, and is typically interpreted as evidence for magnetic disk locking. The eight-day 
bump is not seen for non-disked stars. A two-sided Kolmogorov-Smirnov test indicates a near-zero 
probability (1.277 x 10“^^) that the disked periods vs. non-disked periods were drawn from the 
same distribution. 































































et al. (2001), we identify disk-activity variables via a negative slope between -45 and -90 degrees in 


the H — K,K color-magnitude diagram, and we identify reddening variables with positive slopes 
between 30 and 90 degrees in the same diagramj^ We noted in Rice et al. (2012) that the color 
space trajectories are often linear, and in |Wolk et al.| (2013) we find that many are well-fit using 
this outlier-resistant algorithm. 


4 . 3 . 1 . Examples of typical color behavior 


A few illustrative examples of stars that show well-defined color trajectories are shown in 
Figures [TtI [Tsl and[T^ Figure W shows a fairly typical periodic disked source. The power in each 
band’s periodogram is very high, but the amplitude varies from phase to phase by about 10%. On 
the other hand, the phases themselves seem stable over more than 100 periods as evidenced by the 
points from season two which ht appropriately into the phase-folded light curve. There is limited 
movement in color-color space, but dramatic motion in the direction of reddening (~ 0.5 , ~ 0.7 

Aj ) is visible in the color-magnitude diagrams. The direction of this motion is typically steeper 


than predicted by ISM reddening models (e.g., Rieke &: Lebofsky 1985). 


Figureshows a typical aperiodic reddening object, albeit one with a large amplitude. While 
there is no clear periodic behavior, all the motion in the color-magnitude diagrams still lies generally 
along the reddening vectors. In this case the motion is very close to what one would predict from 
ISM reddening models (e.g., Rieke & Lebofsky 1985). The color-color diagram, on the other 
hand, shows what Rice et al. (2012) called “hybrid motion”, halfway between the expected slope 


of reddening and the CTTS locus. Finally, Figure 19 shows an example of a star which has a 


changing interaction with its disk. This star happens to be a relatively low-amplitude example, 
with AiF ~ 0.2. Its motion in the color-color diagram is almost parallel to the CTTS locus; in the 
K, H — K color-magnitude diagram, it becomes bluer by AH — AT ~ 0.15 as it dims. The J, J — H 
trajectory is best fit with a vector which also gets bluer as it gets fainter, but the scatter is very 
high in this fit and it is hard to determine if there really is a preferred color change in J, J — H. 


4 . 3 . 2 . Aggregate color behavior 


In Figure we show histograms of the htted slopes in color-color and color-magnitude space. 
We include stars of all classes in this aggregate color analysis, and require a star to be at least 
Q = 1 to be included, as well as an additional band-quality criterior The distributions indicate 


In the color-magnitude diagrams, a “positive slope” is down and to the right because the magnitude scale is 
inverted following general astronomical convention. 

^^To be included in the J — H vs. H — K analysis, a star must be Q = 2; to be included in the other color analyses, 
it must be Q = 1 and have high-quality, flag-free data in the bands in question: J and H, or H and K. 
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Fig. 16.— Illustration of two example tracks in color-magnitude space. A “positive” slope (top; 
ONCvar 310 track shown as example) is associated with variable dust extinction, while a “negative” 
slope (bottom; ONCvar 1226 track shown) is associated with accretion variability or changes in 
the inner disk geometry. More variability mechanisms, such as cool and hot spots, exist, and are 
outlined in Wolk et al. (2013). In this figure, the expected direction of dust reddening is derived 


from the extinction law given by Rieke &: Lebofsky (1985); similar results are obtained from more 


modern versions of the extinction law. 
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Fig. 17.— Three-band light curve folded on a period of 2.57 days, plus color-magnitude and 
color-color diagrams, for the variable star ONCvar 500. This star exhibits periodic “AA Tau”-type 
variability (as defined in 16.1) that shows a dust reddening signature in H — K,K space. 
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Time (MJD - 54034.0) H-K 


Fig. 18.— Three-band light curve, plus color-magnitude and color-color diagrams, for the variable 
star ONCvar 539. While the color-magnitude trend is consistent with (and dominated by) red¬ 
dening, the observed deviations from the reddening vector are far in excess of observational errors, 
indicating that more than one process contributes to the measured color changes. 
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Fig. 19.— Three-band light curve, plus color-magnitude and color-color diagrams, for the variable 
star ONCvar 506. This star shows mild color-magnitude variability and color-color variability that 
is consistent with changes in accretion rate or inner disk geometry (Carpenter et al. 2001). 
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a strong bimodal behavior m. K, H — K and weaker bimodality in J, J — H. We interpret these color 
slopes as indicative of disk-activity variability (bluer when fainter) and dust reddening variability 
(redder when fainter). As in [Carpenter et (2001), only stars showing signihcant variability in 
the appropriate bands/color channels that exceeded noise by a factor of 1.5 are shown. The filled-in 
portion of the histogram represents those objects to which we htted color slopes with better than 
20% accuracy. These plots should be directly comparable to Carpenter et al. (2001; their Figure 
20). Our data show a much larger relative fraction of accretion/disk-activity slopes versus dust¬ 
reddening slopes. We find up to 44% of the fitted sources show signs of disk activity, compared to 
between 15-20% by [Carpenter et ^ 2001 


In addition to differences in slope, each variability mechanism has a different maximum vari¬ 
ability amplitude. For example, as summarized in Wolk et al. (2013), while extinction can add 


essentially infinite reddening, cold spots are limited to changes of < 15% and even hot spots are 
limited to changes of about 0.6 mag at K. Disk models are less well-calibrated but appear capable 
of changes of at least 0.75 mag at K. By comparing the observed distributions of color slopes in 
A(J — H)/A{H — K), AK/A{H — K), and AJ/A(J — H), we can infer information about which 
variability mechanisms are most common in our data, as well as identify which color slopes are best 
for distinguishing certain mechanisms. 


4-3.3. Compound variables 

We also note a significant number of stars that show changes in the observed variability mech¬ 
anism. Their trajectories in color-color space and color-magnitude space are not confined to a 
single linear slope, but instead show both “positive-sloped” and “negative-sloped” tracks cleanly 
separated in time, indicating that short-term monitoring may not be sufficient in identifying a given 
star’s primary mode of variability. Indeed, there may be no “primary mode” of variability in these 
objects. We discuss a few examples below: 


ONCvar 494 (Figure [M|), Protostar class, shows a long term rise by over a magnitude in K over 
the hrst season. The rise is not monotonic but has apparent peaks and valleys separated by about 
50 days. The trajectory in color-color space is in reasonable agreement with the reddening vector, 
but not particularly straight. The J, J — H trajectory (not shown) is clearly one of blueing (by 0.6 
in The K, H — K color-magnitude is generally on a blueing trajectory as well, except that 

there are two distinct tracks, one which appeared for the first hundred days or so and a second, 
brighter by 0.2 K, which held for the next ~ 50 days. 

ONCvar 218 (Figure [22]), “No disk” class, displays behavior nearly opposite of ONCvar 494. 
It shows no long-term trend but has significant brightness peaks separated by ~ 50 days. The K, 
H — K color-magnitude diagram is well ht by ISM dust, albeit with a bit of scatter. The color-color 
data is dominated by the very large errors in J band, and does not provide information. 
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ONCvar 212 (Figure [^, Protostar class, was relatively calm for the first ~ 50 days. This is 
followed by an event noticed mostly at J and H as a brightening followed by a significant dimming 
(40% at J), both of which appear reddening-related. The system then brightens gradually at J and 
suddenly at K. The K band event seems to have started between days 130 and 135. The color-color 
track of the large-scale brightening event clearly traces along the CTTS locus. Interestingly, when 
the star is briefly re-observed in season 2 the color-color data are still far to the right of the extended 
CTTS locus, even though the brightness of the star is far fainter in all three bands. 

ONCvar 466 (Figure [2^, Disked class, is an effective counterpoint to the previous figure. Again 
the star’s light curve was relatively flat for the first ~ 70 days. This is followed by an event noticed 
mostly at J and then H but this time only as dimming. As the star reaches its faintest at J and 
H it then brightens in all three bands. The star moves strongly along the CTTS locus starting on 
day 125. At J and H, the star’s behavior is consistent with only a reddening event; the addition 
of K data help flesh out the full picture, where part of the star’s behavior appears to be related to 
disk activity / accretion. 

In addition to well-defined cases of compound color behavior such as those shown above, we 
also find an even larger subsample of stars with significant color variability that shows no clear 
linear trend but could be explained by a complicated mixture of variability mechanisms. 


5. Variability Timescales 


In this section we analyze how various parameters change with the length and number of 
observations in a near-infrared monitoring study. 

We studied how the number of automatically-detected variables (i.e., only variables in quality 
classes Q = 1 and Q = 2] see Appendixj^for details) changed as a function of time baseline (Figure 


At the shortest baselines (~ 5 days) hundreds of variables are detected; this increases to over 800 
after ~ 200 days. When the longest baselines (up to 900 days; not shown in this plot) are included, 
the number reaches the expected total of 857. More than half of the final variables are detected 
within just 5 days of monitoring. 

We studied how the median K amplitude of variables changes as more datapoints are added 
(Figure [2^). At the shortest baselines (< 5 days), the amplitude is below 0.05 mag, but climbs 
steadily to ~ 0.15 mag when the sampling baseline is ~ 200 days. By the end of the ~ 900 day 
study, the median amplitude climbs to ~ 0.16 (not shown in figure). Half of the final variability 
amplitude is reached after about 20 days of monitoring. 


25 i). In the final dataset, there are 857 such “automatic” variables in our catalog of 1203 variables. 


Finally, we studied how the number of “blueing” versus “reddening” stars changed over time 
(Fig. 26). As a class, we refer to these stars as “color variables”. The following plots were made by 
comparing two samples: (a) the original data, un-altered, but sampled over time windows of varying 
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Slope (degrees) 


Fig. 20.— Distributions of the fitted slopes in color-color and color-magnitude space. As in Carpen¬ 


ter et al. (2001), only stars showing rms variability in the appropriate bands or color channels that 


exceeded noise by a factor of 1.5 are shown, with the colored-in bars showing color slopes with bet¬ 
ter than 20% accuracy. Notably, these data show a much larger relative fraction of accretion/disk- 
activity slopes versus dust-reddening slopes than did Carpenter et al. (2001); see Figure 20 in that 
paper. We show in ^ that this disparity is primarily due to our longer observing baseline: the 
longer one observes in the near-infrared, the more accretion/disk-activity variables are seen. 
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Fig. 21.— Three-band light curve, plus color-magnitude and color-color diagrams, for the variable 
protostar ONCvar 494, which shows “compound variability” in color-magnitude space (discussed 
in 14.3). This star exhibits a H — K,K track showing accretion/disk activity from t = 0 to 
t = 80; around t ~ 80, the star gets bluer as it brightens (an indicator that the dust reddening 
has decreased) before resuming the accretion/disk activity until t = 130. At t = 130, the slope in 
H — K,K shifts again, showing variability indicating changes in dust reddening. 
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Fig. 22.— Light curve of the variable star ONCvar 218, which lacks a mid-infrared disk and shows 
“compound variability” in color-magnitude space. This star’s variability moves it back and forth 
along a dust-reddening track in H — K,K space when brighter than K = 15.3, but moves along an 
accretion/disk-activity track when fainter than K = 15.3, especially in the third observing season. 
Its trajectory in the J — H versus H — K color-color diagram is dominated by noise at J and should 
be disregarded. 
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Fig. 23.— Light curves and color-magnitude and color-color diagrams for the protostar ONCvar 
212, which shows “compound variability” in color-magnitude space. This star occupies a “quies¬ 
cent” region in H — K,K space (which is still significantly variable: A{H — K) = 0.2) until t = 130, 
when it abruptly migrates by A{H — K) = 0.6 along an accretion/disk-activity track in H — K,K 
space, brightening by 1.3 magnitudes at K. Between the hrst and second observing seasons, the 
star dims by AK = 0.6, moving along a dust reddening track. In addition to its exceptional color 
behavior, this star shows striking differences between the shape of its light curves at each band. 
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Time (MJD - 54034.0) H-K 


Fig. 24.— Light curve of the variable disked star ONCvar 466, which shows “compound variability” 
in color-magnitude space. Between t = 0 and t = 125, the star moves in H — K,K space along 
a shallow dust-reddening track whose slope may be interpreted as simultaneous reddening and 
disk/accretion activity. Then, at t = 130, the star brightens, moving up in H — K,K space along a 
steeper reddening track. In the third observing season, this alternating behavior is seen to repeat, 
suggesting that this star would show both mechanisms when monitored over even longer timescales. 
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Fig. 25.— (left) Number of automatic variables detected as a function of sampled baseline, (right) 
Median K amplitude as a function of sampled baseline. Error bars depict the median absolute 
deviation (MAD) of variable star amplitudes around the median. 
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length; (b) the data with randomized timestamps, also sampled over time windows of varying length. 
The inclusion of time-randomized data is to distinguish whether the signal simply emerges because 
additional datapoints bring down the noise, or because the addition of data sampled at specific 
times reveals a truly time-varying signal. Multiple non-shuffled samples were created by sliding the 
time window by the size of window considered, allowing us to probe the same temporal space with 
independent subsets of the data. 


As noted in §4.3.2| we find a larger proportion of “blueing” to “reddening” color variables 
in this study compared to Carpenter et al. (2001). By analyzing color trends over progressively 
longer subsets of our time-series data, we find that this difference is due primarily to our longer and 
more complete time coverage (Figure [27]). We found that while 1/3 of the reddening-dominated 
stars could be identified with just a few days of data, at least a month was required to identify a 
similar fraction of disk/accretion dominated objects. A little over one-fourth of the color variables 
never presented a linear slope which could be fitted with precision of better than 20%. Within 
5 days, more than 10 reddening objects were typically found, but less than 5 blueing objects are 
found in these short baselines. By the end of the study, there were 54 “reddening” objects and 
36 “blueing” objects identified; half of the reddening objects were identified within 40 days, but it 
took 90 days to identify half of the blueing objects. There is a substantial difference between the 
shuffled and non-shuffled samples: in the shuffled sample, the reddening and blueing curves have 
a similar shape, and hence a similar time-dependence; but in the non-shuffled sample, the blueing 
curve lags substantially behind the reddening curve and has a lower slope until around 80 days. 
Since the characteristic timescale of blueing objects is more than twice that of reddening objects, 
this is a strong indication that long, contiguous studies are necessary to better understand the 
blueing objects, even as studies around 40 days in length were able to sufficiently characterize the 
reddening objects. 


6. Exceptional Variables 

The tremendous amount of data in a survey such as this has allowed us to identify a number of 
lightcurves of special interest. These systems include “Dippers” and AA Tau variables (^6.1), stars 
which undergo widely-spaced periodic brightenings (^6.2), stars with particularly long periods or 
long-timescale non-periodic variability (^6.3[), and eclipsing binaries (Q. Each is discussed below. 


6.1. Dippers and AA Tau type variables 


Morales-Calderon et al. (2011) note a number of stars that get redder as they get fainter during 


abrupt “dimming” events that are sometimes periodic, sometimes not; they use the term “dippers” 
for these stars. The interpretation of this behavior presented in Morales-Calderon et al. (2011) is 
that these dips are either caused by high-latitude warps in the inner disk (especially in periodic 
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Fig. 26.— Number of reddening (red dots) and blueing (blue X’s) objects as a function of time, for 
the shuffled and non-shuffled samples. There is a very signihcant difference between the shufffed and 
non-shuffled samples, indicating that the different shapes of the non-shufffed red and blue curves 
are due to differences in their intrinsic timescales. Dotted and dashed lines in each color show 50% 
and 100% of the final count of color variables of each type, respectively. Note that the non-shuffled 
blue curve takes much longer to reach its 50% mark than the time-shuffled blue curve. In the 
shuffled data, a few points appear above the “100%” line due to chance selection of non-consecutive 
color points during the shuffling process; this does not affect the interpretation. 

















- 50 - 



Fig. 27.— The number of variables showing significant H — K variability (and a ‘good’ fit to a color 
slope, in the red and blue curves), as a function of observing baseline. By analyzing progressively 
longer subsets of our data, we find that the number of “color variables” increases steadily with 
the length of monitoring. Notably, the disk/accretion variables show the strongest dependence on 
observing baseline: for simulated campaigns spanning less than ~ 80 days, less than half of the 
accretion/disk variables that we find in our final analysis are recovered. The length of our nearly 
continuous first observing season is ~ 180 days; in this plot, the slope of each line becomes shallower 
after this interval, suggesting that if we had nightly observations year-round for the same 894 day 
span, the t > 200 slope of each line would be steeper. Note that the red curve is not monotonically 
increasing because some stars initially show a reddening color trend that is washed out by later 
data (in some cases due to the compound variable behavior outlined in (4.3.3). 
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cases showing AA Tau-type variability) or that the stars are being extincted by overdense clouds 
in the disk atmosphere. In our dataset we also find a number of dippers. 


To clarify terminology, we make a distinction between non-periodic dippers, which may or may 
not show repeated “dipping” events, and the periodic AA Tau-type variables; we use only the term 
“dipper” for the former, and “AA Tau” or “AA Tau-type variable” for the periodic latter class. 
Example dippers are shown in Figures 28 and 29 Example AA Tau-type variables are shown in 


Among individual, non-periodic light curves, a few stand out as deserving particular attention 
as “dippers”. ONCvar 310 (Figure 28) shows a dramatic, well-defined “dip” on top of mostly 
quiescent variability. The quiescent state lasts from t = 0 until t = 135; from t = 135 to f = 160, it 
decreases and then steadily increases in brightness, with A J = 0.92, AH = 0.53, and AK = 0.33. 
The observed wavelength dependence during eclipse closely matches typical ISM exctinction laws. 


Five other prominent examples of one-off or irregular “dippers” are shown in Figure 29 


6.1.1. Periodic “dippers”: AA Tau analogs 


The variable star “AA Tau” exhibits periodic dips in flux consistent with dust-based extinction, 
leading authors to suggest the presence of a magnetically-warped inner dust disk responsible for 
the dips in flux (Bouvier et al. 2003, 2007; Donati et al. 2010). This class of variability appears 
frequently in our data. We identify such variables in our data using three criteria: a) They are 
periodic, as determined by the periodicity criteria presented in 13.2.1 b) Their color variability in 
H—K,K space is significant and consistent with dust-based eclipses c) Their lightcurve morphology 
is not well-behaved (i.e., “cleanly” periodic) enough to be explained by rotational modulation of 
starspots; often, this means large variability amplitudes, an amplitude that changes over time, etc. 
Several example A A Tau type light curves are shown in Figure 30, with a range of periods and 
variability amplitudes. Notably, each AA Tau type variable seems to sweep out a similar narrow line 
(parallel to the CTTS locus of Meyer et al.^1997 ) in J — H,H — K space, even though J — H, H — K 
behavior is not among the selection criteria. 


By applying these criteria to all 539 periodic light curves and investigating them by eye, we 
identify 73 AA Tau variables. This is 6.1% of all variables and 13.5% of periodic variables. The 
class distribution of AA Taus is 2 protostars, 63 disked stars, 3 non-disked stars, and 5 of unknown 
class; clearly this group is by far dominated by disked stars. There are 641 disked variables, and 
229 periodic disked variables in this study. Thus 9.8% of all disked stars are of AA Tau type, and 
27.5% of periodic disked stars are of AA Tau type. 

The period distribution of the AA Tau type stars is shown in Figure Most AA Tau 

type variables have periods less than 15 days, with a concentration below 6 days and another 
concentration between 7 and 10 days. This slight bimodality is reminiscent of the overall disked 
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Fig. 28.— Light curve for “dipper” star ONCvar 310. It shows an abrupt dust eclipse from t = 135 
to t = 160, dimming by 0.9 mag at J and 0.3 mag at K. 


Time 






















- 53 - 


Li gilt curve (MID ^ 54034) 



Q &B tQQ l&B TM KA fiE4 


J-H tfS. H-K 


1.0 




if VS. if-if 
T-r 



AC 


rijUja.4a.'^ai.aiJrt.4 


i.aai.erTi.ui.2iiJ«i.3& 


1.0 

~i—rn— \j^' « 1 

ll.« 

— 1 — 1 — 1 —\— ' 



11^ 


l.L 


lljQ 


Ijfi 




aj 

r ' 

ii.a 

114 

□ 


i].Aeja.4Q.«0Lai.aiJrt.4 
I I I 1.^ j 


a.Bea£iei.a4i.QiLLat.3B 
11.4 



a.>»ja.4a.ia&ai.aijrt.4 



a.EisaiTa.^&Tia.iaa.aE. 
114 


ii.^a 



a.Aaja.4a.iaBLai.aiJt.4 


a.at&a4i.aj iiii^ti.ia 



aja«ja.4a.«0Lai.aiJt.4 



Fig. 29.— Five more “dipper” light curves, each somewhat irregular. Columns from left to right: 
K band light curve; J — H vs. H — K color-color diagram; K vs. H — K color-mag diagram. 
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period distribution from the middle panel of FigureWhen plotted against median K magnitude 
of each star (Figure|3^), we see a strong trend for fainter stars to have shorter periods. The scarcity 
of A A Taus in the bottom-right quadrant of this plot is primarily a selection effect, explained by 
our low sensitivity to variables at these periods and brightnesses: A noisy signal is easier to confirm 
observationally if it repeats many times. Since bright stars have less observational noise than faint 
stars, their periods are more easily detectable even when the periods are long. The worse signal- 
to-noise for fainter stars means more repetitions are needed to confirm a signal. The lack of AA 
Taus in the top-left quadrant, however, is not explained by a selection effect, and may indicate 
that AA Taus of a given period have a maximum luminosity. This is perhaps explained by noting 
that the brightest stars will sublimate the dust that would orbit at small distances, thus putting a 
minimum period for bright stars. As discussed in Section the physical properties of such young 
stars and their disks have been measured by Eisner et al.| ( |2009 ): the inner disk walls are in the 
range 0.14 — 0.18 AU, with disk wall temperatures of 1260 — 1850 K. 


6.2. Long-Period Brighteners 


Most periodic variables either a) spend equal amounts of time brighter and fainter than the 
mean (e.g., starspots or any nearly sinusoidal variability) or b) spend most of the time brighter 
than the mean (e.g., eclipsing binaries or stars that undergo occultation). In addition to these 
types of variables, we observe three periodic variables with cleanly “inverted” light curves, i.e., 
they spend most of their time fainter than their mean brightness; four more stars show similar 
behavior, but not as cleanly, and are discussed in Section |6.3[ Among YSOs, this general class of 
near-infrared behavior has been previously seen in the long-period variables CHS 7797 in the Orion 
Nebula Cluster (17.8 day period; Rodrfguez-Ledesma et al.|2012[ 2013), WL 4 in p Oph (130.87 day 
period; Plavchan et al.||2008 ), and YLW 16A also in p Oph (92.6 day period; Plavchan et al.||2013 ), 
with KH 15D in NGC 2264 (48 day period; Kearns &: Herbst||1998 Windemuth Sz Herbst||2014 and 
references therein) as an optical counterpart. Importantly, all such variables show periods greater 
than 15 days, suggesting dynamics longer than either typical stellar rotation or the orbital period 
of the inner disk. In this dataset, we find three unambiguous examples of long-period brighteners: 
the previously noted CHS 7797, known here as ONCvar 497 (Figure]^, ONCvar 149 (36.4 day 
period; Figure [33|), and ONCvar 1226 (88.5 day period; Figure [34]). This type of periodic behavior, 
sometimes called “inverse eclipses” (e.g., by Parks et al. 112014]), has been suggested to be due to 


occultations of binary stars by their circumbinary disks (Parks et al. 2014, Plavchan et al. 2008 


2013). Such systems may provide key insights into circumbinary disks, where circumbinary planets 


such as Kepler-16b (Doyle et al. 2011) would be expected to form. 


KH 15D in NGC 2264 has a period of 48.37 days and is interpreted as a pair of T Tauri stars 


surrounded by an inclined, warped, disk-like ring of non-ISM dust ( 

Windemuth &: Herbst 

2014) 

Its periodic behavior was discovered in 1998 ( 

Kearns & Herbst 

1998 

) and it has since been well- 

characterized at many wavelengths, with polarization and spectra. 

Windemuth &: Herbst 

(2014 
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Fig. 30.— Ten example AA Tau light curves 
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note the possibility of a young planet shepherding the well-defined ring structure and possibly 
contributing to J and H excess flux seen toward the system. 


These long-period brightening sources are reminiscent of peculiar stars seen in the embedded 
star forming region Lynds 1688 in p Oph. WL 4 and YLW 16A in p Oph, the 131-day and 93-day 
periodic systems found, respectively, by Plavchan et al. (2008) and Plavchan et al. (2013), and 
discussed further by Parks et al. (2014), have near-infrared behavior remarkably similar to the 
examples displayed below (especially ONCvar 149 and 1226), so the analysis presented in these 
papers may serve as a useful guide for interpreting the two ONCvar objects. In summary, the 
lightcurves and SEDs of these stars are well-modelled by a triple system including an inner binary 
pair eclipsed by their circumbinary disk; both systems have resolved the outer tertiary companion. 
The color behavior during the “bright” part of the lightcurve gets redder at J — H and H — Kg in 
each system. 


6.2.1. ONCvar 4 79 (CHS 7797) 


ONCvar 479 (CHS 7797; Carpenter et al.|2001 Rodrfguez-Ledesma et al.|2012 2013) shows 


an 


unusual periodic brightening cycle with period ~ 18 days and a large amplitude of ~ 1.7 mag. This 
star has been characterized extensively for photometric variability in optical and infrared ([3.6] and 
[4.5] Spitzer bands) as well as spectroscopically in two phase points in its periodic light curve. This 
star’s period was accurately derived to be 17.786 days by |Rodrfguez-Ledesma et al.| (2012) using 
six years of /-band monitoring; its variability was also characterized at R, z', [3.6] and [4.5] bands, 
showing variability amplitudes ~ 1.7 in /?, / and z' and decreasing slightly at longer wavelengths. 
The star’s SED indicates that it is a Class I object, and its spectral features show that it is a 


low-mass object near the stellar/substellar boundary (Rodrfguez-Ledesma et al. 2013). The star 
was noted to show colorless variability at wavelengths A < 2pm, but at A > 2pm it became redder 
when fainter; this was interpreted as evidence for dust grains that had grown to 1 — 2 pm. In 


Rodrfguez-Ledesma et al. (2013), the authors identify spectral differences between two light curve 


phases that suggest the system has a similar geometry to KH 15D. 

In our observations (Figure [3^, the star varies with amplitude A/L ~ 1.5, AJ ~ A// ~ 1.4, 
within the fluctuations seen in the optical monitoring, and its flux changes are on average colorless. 


which is consistent with the Rodrfguez-Ledesma et al. (2012) study. The period we derive from 


these observations is 17.795 days, consistent within 0.01 day of the six-year monitoring study. 
Many of our JHK observations are contemporaneous with the six-year monitoring of this star at 


/ (Rodrfguez-Ledesma et al. 2012) 
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6.2.2. ONCvar U9 


Another long-period brightener is ONCvar 149 (Figure [33]), which is redder and has a much 
longer period of 36.40 days. Its Megeath et al. ( 2012[ ) classification is “protostar” and it is highly 
reddened, with H — K varying around 2.3. Due to this star’s faintness at shorter wavelengths, our 
data contain only one measurement at J (J = 19.90 ± 0.24). Even though this star’s faintness at J 
and H means that its measurements at those bands have high uncertainties, it is still clear that the 
variability amplitude is much higher at K (one cycle: /S.K = 0.64) than at H (one cycle: AH < 0.3), 
in stark contrast to the colorless or weaker-with-long-wavelength variability seen in ONCvar 479 
(CHS 7797). This K, H — K slope indicates that disk/accretion activity dominates the variability. 
The differing wavelength dependence of variability may indicate that the superficially similar light 
curve morphology seen in ONCvar 479 (CHS 7797) and ONCvar 149 does not ultimately originate 
from analogous physical conditions. 


6.2.3. ONCvar 1226 


ONCvar 1226 is another variable showing “long-period brightener” light curve morphology 


(see Figure 34). It shows the characteristic periodic bursts that fold on a steady period with 
modulations in amplitude. This light curve is particularly notable for two reasons: a), its very 
long 88.5 day period is longer than previous NIR monitoring campaigns of the ONC; and b), its 
bursts get stronger with increasing wavelength ~ they are clearly seen at K, and show an even 
larger amplitude at [3.6] and [4.5] bands, but the J band light curve has weak anticorrelated (i.e., 
fading) behavior instead; this anticorrelation gives ONCvar 1226 a Stetson index of 0.97, just below 
the nominal 1.0 cutoff. This also sets ONCvar 1226 clearly apart from ONCvar 479 (CHS 7797), 
whose light curve in bands shorter than 2;um, have similar shape and amplitude including R, I 
and z'. Unlike ONCvar 149, which is faint and highly reddened [K ~ 15, J ~ 20), making optical 
followup impossible, ONCvar 1226 is bright and relatively unreddened {K ~ 11, J ~ 12), so it 
would be an excellent target for optical spectroscopic followup, as was done with CHS 7797 in 
Rodriguez-Ledesma et al. (2013). 

Based on the 88.5 day period, ONCvar 1226 is expected to have undergone ten brightening 
events during our 894-day observing window. Five of these overlap with our actual observing 
epochs, and all five are confirmed in the light curve. (This is fortuitous, as our 120 nights have 
only ~ 13% coverage of the 894-day window, and yet we catch 50% of the brightenings.) The K 
amplitude of each brightening event varies considerably, from AK ~ 0.11 to 0.3, but most are near 
0.3. This star has not been previously noted to be periodic or exceptional, but fortunately, the 


public Spitzer YSOVAR light curve data in the Orion Nebula Cluster taken in late 2009 (Morales- 


[Calderon et al. 2011) coincide with the third brightening event after the last one we observed 
(Figure [ 3 ^. We accessed the YSOVAR data after finding the star to be of interest from our data, 
and were able to further confirm the ~ 88.5 day period using these IRAC data, which show the 
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latter half of the star’s brightening event down to its quiescent level in the [3.6] and [4.5] IRAC 
bands, with amplitude 0.38 and 0.45 respectively. 

At K band, hints of possible “secondary” brightenings at half-phase are seen around t = 167 
and t = 876 with amplitude AK ~ 0.08, but the data are not conclusive. A corresponding half¬ 
phase brightening is not seen at f = 79, where one might be expected. Notably, the IRAC [3.6] and 
[4.5] data from YSOVAR each indicate a clear, plateau-shaped rise in brightness of order 0.06 mag 
(well above the photometric noise limit) at half-phase (Figure [35]). 

Our JHK data could also be consistent with a period twice as long, or 177.0 days; the three 
cleanest brightenings (at t = 35, t = 389, and t = 743 all fall on each other when folded in this way. 


6.3. Long Periods and Long-Term Variability 
6.3.1. Periodic 


In addition to the long periods seen in [6.2 (ONCvar 149 and 1226, Figures 33 and 34), 12 


other stars with near-IR trends indicating cycles longer than ~ 30 days are apparent in our data 
(see Figure [ 3 ^ for several examples). The aforementioned stars ONCvar 149 and 1226 have periods 
likely to originate in the orbital period of a binary companion; in other cases, we speculate that 


slow-rotating spotted stars induce a clean long-term period (such as in Figure 14), while in other 


cases a more messy period is induced by dynamic behavior associated with the inner edge of a disk 
interacting with its star (this is suggested because such systems often show “blueing” behavior in 


K, H — K color-mag space). The observations presented in Plavchan et al. (2008) and Plavchan 


et al. (2013) represent some of the first identifications of long-period systems. 


We restricted our automated period search in §3.2| to periods below 50 days in order to guaran¬ 
tee that the observed periods were reliably repeatable, as many irregular variabilities in young stars 
may appear sinusoidal over some time frame only to change its variability pattern later. In this 
section we focus on targets whose apparent periodicity is nonetheless worthy of comment, although 
we cannot constrain the absolute frequency of such systems. 

Some stars show what might be classified as purely irregular variability that, observed over 


long term, shows an overall periodicity. ONCvar 637 (Figure 36, row 1) is a system showing 


quasi-periodic behavior over a 42.5 day period, with varying amplitude each cycle, in which the 
variability is dominated by changes in K excess. Such a system is inconsistent with A A Tau- 
type dust-occultation behavior, and the 42.5 day period is also longer than the expected rotation 
period of the inner disk that emits at K band. The period may be caused by the presence of a 
binary companion, perhaps modulating the accretion rate or the disk structure via gravitational 
interaction with an orbiting point mass. Note that ONCvar 637 avoids a definitive “long-period 
brightener” classification as its “bright” state is not narrow compared to its “faint” state, and also 
its light curve is somewhat irregular, but its period and color behavior are otherwise comparable to 














- 59 - 



AA Tau period (days) AA Tau period (days) 


Fig. 31.— (left) Distribution of periods among AA Tau type variables, (right) Median K magnitude 
of AA Tau type variables versus period. The apparent trend in the data is likely an observational 
bias, as long periods are intriniscally easier to detect for bright stars. 
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Fig. 32.— Light curve for ONCvar 479 (also known as CHS 7797), an example of a “Long-Period 
Brightener”. 
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Fig. 33.— 


Light curve for ONCvar 149, an example of a “Long-Period Brightener”. 
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Fig. 34.— Light curve for ONCvar 1226, the longest-period example of a “Long-Period Brightener” 
found in this dataset. 
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Time {MJD ^ 54034) 

Fig. 35.— 40-day YSOVAR light curve for ONCvar 1226, showing a brightening event and possible 
secondary brightening roughly 7 months after the end of our UKIRT monitoring program. The 
[3.6] and [4.5] data show a continuation of the 88.5 day period seen at H and K. 



Table 6. Long period periodic stars 


ONCvar ID 

Period (d) 

LPB 

231 

70.48 


479 

36.40 

Y 

529 

92.51 

C 

637 

42.5 

C 

665 

44.27 

C 

801 

96.21 


802 

158.8 

C 

816 

92.35 


819 

37.0 

C 

826 

30.22 


843 

75.16 

C 

874 

72.88 

C 

1027 

37.82 


1226 

88.5 

Y 


Note. — A “C” in the last 
column denotes stars that are 
identified as candidate long-period 
brighteners, while a “Y” denotes a 
definite identification as such. 
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the long-period brightners presented in the previous section (see Figure 34 and 33). ONCvar 874 
(Figure [36| row 3) shows similar behavior during a roughly 50 day period, in which the amplitude 
modulates over the course of our observations, but the color signature and qualitative morphology 
is the same as in ONCvar 637. ONCvar 665 (Figure]^ row 4) is a third example of a star showing 
this approximate behavior, in which irregular variability dominated by an accretion/disk track in 
H — K,K space is modulated by a period around 40-50 days. This one has a 44.27 day period. 

Some stars vary smoothly over such long terms that not even our long time baseline can 


confirm that they are repeatably periodic. ONCvar 663 (Fig. 37) does not show definitively 
periodic behavior, but shows nearly sinusoidal long-term changes continuously throughout our first 


observing season. The light curve of ONCvar 802 (Fig. 36, row 2), in contrast, has a believable 
158.8 day period, where the modulation seen over the course of the hrst season begins to repeat 
and the second-season observations fall on top of season one when folded on 158.8 days. 


6.3.2. Non-periodic 

25 stars show dramatic long-term variability that is not periodic. 

ONCvar 957 (Figure [38]) shows a clean period at J but a long-term trend at K. The most 
natural way to explain such strange behavior is to assume that the J-band short-period behavior 
originates on the stellar photosphere in this system, analogous to what would be observed at optical 
wavelengths, while the 77-band long term variability originates in the disk. Indeed, the K variability 
seen in ONCvar 957 shows an unambiguous disk/accretion track in H — K,K space. 

Further shown in Figureare ten examples of prominent, high-amplitude long-term variables 
which do not show strong signatures of periodicity. 

Finally, it is worth noting that many or most of the non-variable objects monitored in this 
survey could show non-periodic variability on timescales longer than 900 days, if only they could 
be monitored for long enough. 


7. Eclipsing Binaries 


Eclipsing binaries are astrophysically important objects because they can yield fundamental 
stellar parameters such as masses, radii, and temperature ratios. Our survey is suited to find 
eclipsing binaries in the direction towards the ONC. In the field, less than one percent of stars are 
eclipsing binaries, but young stars have considerably larger radii than on the main sequence, and 
hence the chance that PMS binaries show eclipses should be correspondingly higher. Our survey 
has monitored several thousand stars, so we would expect to find a number of eclipsing binaries. 
Indeed, we uncovered 10 eclipsing binaries in the area studied, six of which were previously known 


(shown in Figure 40). This is in fact fewer than we had expected given the larger radii of PMS 
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stars. In principle, the paucity of eclipsing binaries might result if the in-spiral phase required to 
form a close binary would on average take longer than the mean age of the current population of 
the ONC. Alternatively, this could, at least in part, be a detection problem, since the irregular 
variability of many of the stars in our study could impact our ability to recognize all but the most 
pronounced eclipses. An even longer monitoring campaign may be required to identify additional 
eclipsing binaries in the ONC. 


In the following we discuss the four newly discovered systems. Unfortunately, we have been 
unable to determine conclusively whether the four new eclipsing binaries are young. Two (V1916 


Ori and V2716 Ori) were reported to show low-level irregular variability by Carpenter et al. (2001); 
marginal levels of variability were seen out-of-eclipse in our data, not exceeding about 0.1 mag in 
each band. None of the four systems show evidence of strong infrared excess based on near-infrared 
and WISE photometry. However, it is perhaps not surprising that a young eclipsing binary would 
not have much infrared excess, depending on the prior orbital evolution after the binary was formed, 
and many of the known PMS eclipsing binaries similarly do not have infrared excess, e.g., JW 276, 
JW 380, VI174 Ori, and Par 1802. Spectroscopy may be required to settle whether these four 
systems are young or only line-of-sight associations with the ONC. 

In Table we list the four newly discovered eclipsing binaries and their basic data, together 
with nine other eclipsing binaries that fall within our WFCAM field. Of these nine, three are 
saturated in all three filters in our data, but for the other six our data conhrm them as eclipsing 
binaries. Eight of these eclipsing binaries are proven young objects (e.g., Stassun et al.||2004 2006 


Irwin et al. 2007; Cargile et al. 2008, Morales-Calderon et al. 2012). 


The period fits for these eclipsing binaries were obtained using the Fast Chi-Squared method 
(Fy^) discussed in Section 3.2 Unlike the Lomb-Scargle periodogram, which searches for sinusoidal 


signals and is not sensitive to periods below twice the sampling rate (here, 2 days), the Fx^ method 
is able to recover periods for non-sinusoidal light curves even below two days. 


7.1. V1916 Ori 


This new eclipsing binary has a short period for which our best estimate from the limited 
observations is 1.112403 days, seen in Figure \4l\ top row. The secondary minimum is at phase 
0.5, indicating a circular orbit. The minima appear to be slightly unequal, suggesting a surface 
temperature difference of a few hundred degrees. A variable star name was given to this star, since 


Carpenter et al. (2001) found the star to be consistently variable throughout their observations. 


with amplitudes in the range 0.11-0.12 and a Stetson index of 1.34. The variability we see is much 
less than that, and only marginally larger than expected from noise. If real, perhaps the variability 
seen by Carpenter et al. was due to transient spots. The 2MASS colors suggest little or no infrared 
excess. The star was also observed by WISE, and when placed in a Wl—W2 vs W2—W3 color-color 
diagram it falls in an area where very few young stars are found (Koenig et al.|20T2). Further study 
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is required to establish whether the star is young or not. 


7.2. ONCvar 1147 


This newly discovered object has a characteristic W UMa-type light curve (see e.g., |Lucy|19^ 
with no constant parts, indicating that the stars are highly deformed. The period is extremely short, 
only 9.492720 hours, which accounts for the deformation of the stars. The light curve has been 


modelled and Figure 42 shows a model ht to the J-band light curve. The solution indicates a 
circular orbit with an inclination of 74.3° and a mass ratio of unity. The sum of the radii of 
the stars is 0.75 times the semimajor axis, indicating that the system is very near contact. The 
secondary minima get gradually deeper relative to the primary minima from J over H to K. The 
near-infrared colors show no evidence for an infrared excess and little or no extinction. The star 
was also observed by WISE, and when placed in a Wl—W2 vs W2—W3 color-color diagram it falls 


in a region where few young stars are found (Koenig et al. 2012). The star’s lightcurve is shown 
with color information in Figure [dTj middle 


row. 


7.3. ONCvar 1190 

This newly discovered eclipsing binary has poor phase coverage, with only half a dozen ob¬ 
servations within the eclipses, as seen in Figure \4^ bottom row. Nonetheless the coverage of the 
constant parts of the light curve is good, and the data show well-separated narrow minima indicat¬ 
ing non-interacting stars. Our best estimate for the period is 3.031 days. The star is not included 
in the 2MASS point source catalog, although it is visible on individual frames. The position of the 
star in the J — H vs. H — K color-color diagram indicates no evidence for circumstellar material, 
only a considerable extinction, in fact the star is not optically visible on a SERC /-band sky survey 
image. We can thus not exclude that this is a background held star unrelated to the ONC, although 
it is located close to young stars, the nearest young star (V1775 Ori) is only 23 arcsec away. 


7.4. V2716 Ori 


While the eclipses of this object are pronounced with narrow deep minima up to half a mag¬ 
nitude, the eclipses are poorly covered; see Figure [4^ The minimum spacing between the observed 
eclipses is about 20 days, and the period of 3.22 days is highly uncertain. The star has a GCVS 
variable star name because Carpenter et al. (2001| found long-term variability during their survey. 


with a Stetson index of 1.18. We see little evidence for variability outside eclipses. The infrared 
colors give no evidence for infrared excess, and the system suffers little or no extinction. 
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Fig. 36.— Composite light curve of four long-period variables that we identify as candidate long- 
period brighteners. Their color behavior is quite consistent, and similar to the color behavior of 


the dehnite long-period brighteners shown in Figures 34 and 33 


Table 7. Eclipsing Binaries towards M42 within WFCAM held 


ONCvar 

Other ID 

2MASS 

0:2000 

S 20 QO 

1“ 

J“ 

H“ 

K“ 

P[d]*’ 

Ref.= 

47 

V1916 Ori 

J05340205-0517256 

05:34:02.06 

-05:17:25.7 

17.30 

15.29 

14.18 

13.72 

1.11 

this paper 

1147 


J05362022-0514214 

05:36:20.21 

-05:14:21.7 

18.18 

16.74 

16.19 

15.97 

0.40 

this paper 

1190 


J05365393-0525550 

05:36:53.94 

-05:25:55.1 


16.67 

15.24 

14.61 

3.03 

this paper 

1191 

V2716 Ori 

J05365496-0458198 

05:36:54.97 

-04:58:19.9 

16.18 

14.27 

13.73 

13.54 

3.22: 

this paper 

122 

V1174 Ori 

J05342785-0541378 

05:34:27.85 

-05:41:37.8 

12.36 

11.16 

10.57 

10.39 

2.63 

Sta04 

427 

JW 276 

J05350571-0523540 

05:35:05.71 

-05:23:54.1 

13.66 

12.10 

11.47 

11.15 

20.48 

Morl2 

518 

Par 1802 

J05351113-0536511 

05:35:11.14 

-05:36:51.2 

12.44 

11.04 

10.24 

9.94 

4.67 

Car08 

1227 

JW 380 

J05351214-0531388 

05:35:12.15 

-05:31:38.9 

13.55 

12.06 

11.44 

11.15 

2.65 

Irw07 


V1501 Ori 

J05351554-0525140 

05:35:15.55 

-05:25:14.1 

11.55 

10.02 

9.05 

8.73 

? 

Morl2 


Ori E 

J05351577-0523100 

05:35:15.76 

-05:23:09.9 

6.63 

6.64 

6.24 

6.06 

9.89 

Morl2 


[H97b]9209 

J05351802-0522052 

05:35:18.03 

-05:22:05.4 

14.69 

11.69 

10.43 

9.54 

5.62 

Morl2 

788 

V2384 Ori 

J05352184-0546085 

05:35:21.84 

-05:46:08.5 

17.28 

14.64 

14.02 

13.53 

9.78 

Sta06 

1100 


J05360595-0500413 

05:36:05.95 

-05:00:41.2 

17.60 

14.90 

14.11 

13.80 

3.57 

Morl2 


^/-magnitudes are from the literature, J, //, A'-magnitudes are mean values from our photometry, or, in the case of the three stars 
V1501 Ori, 0^ Ori E, and [H97b]9209, which are all saturated in our data, from the 2MASS point source catalog or |Carpenter et al.| < |200l] | 

Periods rounded to two decimals. Periods for the four new eclipsing binaries are best estimates from the current data; a trailing colon 
on (3.22:) indicates a tentative period fit. 

‘^CarOS: Cargile et al. (2008), Irw07: Irwin et al. (2007), Morl2: Morales-Calderon et al. (2012), Sta04: Stassun et al. (2004), Sta06: 
Stassun et al. (2006), 
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Fig. 37.— Light curve for slowly-varying star ONCvar 663. It seems to show continuous cyclic 
variations, but our observations are unable to constrain whether its lightcurve repeats regularly. 
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Fig. 38.— Light curve of ONCvar 957, which shows different qualitative behavior at different 
wavelengths; at J, it is neatly periodic, but at K, it has a long-term rise in brightness perhaps 
associated with a steady increase in accretion rate. 









































68 


11.7 
114 
114 
174 
17.1 

17.7 
17.7 
17.4 


Ught curve (MJD - 54034) 




- 

rV«r^ 


I i^ONCvar ^ 


'w‘ 

I' 1 



114 

174 

17.7 

174 

17.4 

174 

174 


100 1S4 701 7S4 0M 


■ j, ' '' 


I > ONCvar1156 


I. I 


100 ISO 701 7S4 404 



144 

14.7 

144 

144 

144 

1S4 

15.7 




I I ONCvar 605 

V 

W II 


' • -s 




A' VS. If-K 


—.■Vi 

& 


04040.40 40i81.0141.4 

I I I |j 


0.00 l.M 1.141.771 441.60 

4 * a _ 



100 ISO 701 7S4 


04040.40 40181.01.71.4 


0.700760470i480S40.60 


Fig. 39.— Ten light curves of long-term non-periodic stars (labelled in panels). Columns from 
left to right: K band light curve; J — H vs. H — K color-color diagram; K vs. H — K color-mag 
diagram. The stars are roughly organized hy H — K, K color-magnitude slope: stars that redden 
when faint are sorted to the top, followed by stars that become blue when faint. 
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Peroti-folded light curve 
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Fig. 40.— Six previously-discovered eclipsing binaries that have light curves in our data. 
Table for period references. 



Fig. 41.— Three of the four newly discovered eclipsing binaries. 
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Fig. 42.-— The newly discovered eclipsing variable ONCvar 1147 shows a W UMa-type light curve, 


indicating highly deformed components near contact. A WD model (Wilson & Devinney 1971 


Wilson 1979, 1993; Vaz et al. 2007 and references given therein) of the system is here laid over 


the J-fTAT-light curves. The K magnitude values shown were displaced by —0.3 mag, in order to 
prevent overlap with the if-light curve. 
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Fig. 43.— One of the four newly discovered eclipsing binaries, with a tentative period of 3.223 
days. 
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8. Concluding Commentary 


In addition to the previous studies of infrared variability in the Orion Nebula Cluster carried out 


by Carpenter et al. (2001) and Morales-Calderon et al. (2011), a number of other recent variability 


studies of young stars have been published. In this section we relate our results to other such 


studies, which include near-infrared monitoring in p Oph (Parks et al. 2014), and Cyg OB7 (Rice 


eTalllM^ [Wolk et al.||2013D , as well as Spitzer mid-infrared monitoring in NGC 2264 (|Cody et al. 


2014; Stauffer et al. 2014) and IC 348 (Flaherty et al. 2013), Additionally, disk-related long-term 


variability in young stars was studied by Findeisen et al. (2013) in the North America Nebula at 
R band. 

Carpenter et al. (2001; JHKg) and Morales-Calderon et al. (2011; [3.6] and [4.5]) made the first 


attempts to get a detailed handle on variability in their respective wavelength regimes. Carpenter 


et al. (2001) monitored ~ 18,000 stars at JHKg in a wide region overlapping the current study, 


and found 1235 variables, noting that much of the variability seen could be explained by cool or 
hot starspots, with some contribution from variable dust obscuration and accretion/disk activity. 


Morales-Calderon et al. (2011) monitored 2069 objects (1249 known Orion YSOs plus 820 other 
likely Orion members) and found 1259 variables; they took special note of a class of stars which 
showed periodic or aperiodic flux dips consistent with short dust eclipses, called “dippers”. [Flaherty] 


et al. (2013) found 125 variables among 269 known cluster members in IC 348, and used models 


of the inner disk structure to conclude that at mid-infrared wavelengths, structural perturbations 
(kinks, warps, and dents) in the inner disk could drive much variability over periods of days to weeks. 
In NGC 2264, Cody et al. (2014) and Stauffer et al. (2014) monitored ~ 1500 stars including 490 
cluster members and > 1000 field stars; Stauffer et al. (2014) found 23 “burst-dominated” light 


curves among them, while Cody et al. (2014) studied the variability of 162 T Tauri stars, presenting 
“stochasticity” (or departures from perfect periodicity) as an important metric when classifying 


YSO light curve morphology. In the North America Nebula, Findeisen et al. (2013) monitored 


~ 59,000 stars and found 3000 variables, including 41 distinct “bursters” or “faders”. [Rice et ^ 
(2012) studied ~ 9200 stars towards Cyg OB7 and identified ~ 200 variables, of which 30 showed 


a near-infrared excess; using the Cyg OB7 data, Wolk et al. (2013) built on the work of Carpenter 


et al. (2001) to identify physical mechanisms that could be associated with the observed color 


trends. 

An advantage to our focus on a region such as the Orion Nebula Cluster, with its rich history 
of detailed studies, is that we have independent determination of the evolutionary classes of nearly 


all of our variables (primarily via Megeath et al. 2012). Thus we can comment directly on the 


near-infrared variability characteristics of young stars in different classes. Morales-Calderon et al. 


(2011) monitored a mostly overlapping region and also used the same data for evolutionary class 


assignment, so we can directly compare our near-infrared variability results to their mid-IR results, 
with the caveats that our study has a much longer time baseline. Of their 126 Class I sources, 
106 were variable, 71 with amplitude > 0.2 mag (67%), 15 with amplitude > 0.5 mag (14%), 
3 with amplitude > 1 mag (3%), and 3 with [4.5] — [3.6] color amplitude > 0.2 (3%). Among 
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Q = 2 stars in this study, we identify 11 Class I stars, all of which are variable: 9 with amplitude 
> 0.2 mag (82 ± 27%), 6 with amplitude > 0.5 mag (55 ± 22%), 3 with amplitude > 1 mag 
(27 ± 16%), and 8 with H — K color amplitude > 0.2 (73 ± 26%). These numbers suggest two 
different interpretations: perhaps the most extremely variable Class I stars are also brighter in the 
near-infrared than less-variable ones (reflecting an observational bias in our data), or that their 
measured variability increases dramatically over longer timescales. 


Carpenter et al. (2001) measured the fraction of significant color variables that showed either 


a reddening behavior (associated with variable dust extinction) or blueing behavior (associated 
with accretion/disk activity); they found ~ 6 times more reddening-type stars than blueing-type 
stars. In this study we identify a ratio of 1.5 times as many reddening to blueing stars ( ^4.3.2 ), 
thus finding substantially more blueing-type relative to the previous study. Part of this difference 
is likely to be from the fact that the region we study is smaller, more focused on the youngest stars 
in Orion, and deeper at H and K by several magnitudes, giving us more sensitivity to the types of 
stars likely to show large color changes. But a large part of it comes from the fact that our observed 
baseline is longer: we analyze in the number of color variables detected as a function of time, 
and find that while the number of reddening sources rise quickly to their hnal value, the number 
of blueing sources takes longer to rise, indicating a longer intrinsic timescale for the underlying 
variability mechanism. This shows that monitoring programs longer than 1-2 months are necessary 


to accurately study variability, and that the conclusions reached by Carpenter et al. (2001) on the 


nature of infrared variability in young stars are only true on shorter timescales. This may reveal 
a population of disks whose structures are quasi-unstable and that these systems take 1-2 months 
to reveal themselves. A possible origin of such instability may lie near, but not at, the inner disk 
boundary. 


While infrared monitoring studies such as Carpenter et al. (2001), Morales-Calderon et al. 


(2011), Flaherty et al.| (2013), and Cody et al. (2014) are limited by their relatively short contiguous 
time coverage (~ 1 — 2 months), some longer continuous or nearly-contiguous infrared surveys of 
young stars have been carried out, especially Parks et al. (2014; 2.5 years in p Oph), Rice et al. 
(2012; 1.5 years in Cyg OB7), and the current study (2.5 years). These longer studies have enabled 
a better understanding of changes in young stars on timescales significantly longer than the stellar 
rotation period or the orbital period of the inner disk; since infrared studies are best suited to 


studying disk activity, this longer timescale is well warranted. Rice et al. (2012) showed that over 


1.5 years, 36% of stars that show a iF-band excess will drift in color-color space between “having 
an excess” and “not having an excess”. Parks et al. (2014) show that some stars have aperiodic 


variability with timescales between 64 and 790 days. They also discuss seven stars in p Oph which 
clearly demonstrate multiple variability mechanisms over their ~ 2.5 year data. They focus on stars 
which show different, overlapping periodic signatures with different amplitudes and/or morphology 
(such as rotation versus orbital disk phenomena), or periodic versus aperiodic variability in the 
same light curve 


Similar behavior is seen in this study, e.g., ONCvar 957 ( ^6.3[ Figure 38); 
in addition, we note individual stars which show different variability mechanisms that impart 
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distinct, orthogonal color behavior in H — K, K space at different times ( ^4.3.3 ). In this study we 
also identify long-period systems, and demonstrate that the timescales involved for disk/accretion 
variability is typically longer than many previous YSO infrared monitoring studies. In addition 


to these infrared studies, the optical study done by Findeisen et al. (2013) lasted 3.3 years and 
focused on disk-related variability. They study a number of “bursters” and “faders”, finding that 
fading behavior is twice as common as bursting behavior, though individual objects change the 
duty cycle of their bursts and fades from year to year. We now know that JHK variability, 
unlike most optical variability, happens on many different timescales: while “early” studies focused 
on variability tied to the rotation of the star and the orbit of the inner disk edge, these more 
recent results, especially from the present study, indicate that significant changes in accretion rate 
and disk structure occur on one- to many-month timescales. If much of this month(s)-timescale 
“blueing” variability is primarily accretion-based, it may imply that a significant fraction of disks 
have density enhancements in the inner 0.1 — 0.2 AU capable of influencing accretion. Such density 
enhancements have typically not been considered in previous models of young stellar variability, 
which usually consider stable accretion. Overall, these examples of “non-stationary” variability 
( jFeigelson &: Babu 2012) indicate that short-duration monitoring is not sufficient to characterize 
the near-infrared variability of individual pre-main sequence stars, as many stars change the typical 
amplitude, morphology, and mechanism of their variability over timescales of months to years. 

The conclusions of this analysis of near-infrared variability in the Orion Nebula Cluster over 
cl 900 day span: 


1. We have identified 1203 variable stars in a nearly square degree in the Orion Nebula Cluster; 
among an unbiased subset of these, about 24% are periodic variables. 

2. Guided by infrared and X-ray evolutionary class determinations from Megeath et al. (2012) 
and Pillitteri et al. (2013), we find that NIR variability amplitude in all bands and colors is 
a strong function of class, with most variability detected in the least evolved stages. 

3. Color behavior in H — K,K color-magnitude space contains information about physical pro¬ 
cesses behind variability, with the dominant color-variable processes involving dust extinction 
and accretion/disk activity: 

• We observe 1.5x more color variability among YSOs due to dust extinction than accre¬ 
tion/disk activity, but the proportion of color variables showing accretion/disk activity 
increases steadily with observed time baseline, as does the total number of color vari¬ 
ables. 

• A substantial number of stars show “compound” color variability in which multiple 
modes of color behavior are seen separated in time for an individual star. This suggests 
that many more stars may undergo such behavior over longer timescales, and hints that 
perhaps “primary mode of variability” is appropriate only for shorter time-intervals in 


some cases. 
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4. We identify 73 “AA Tau” type periodic variables, which constitute 6.1% of variables and 
13.5% of all periodic variables. They are predominantly disked stars; 9.8% of all disked stars 
exhibit periodic, AA Tau type variability. 

5. We identify a number of long-term periodic variables, including three “long-period bright- 
eners” in the class similar to the 130-day and 93-day periodic objects identified in Plavchan 
et al. (2008, 2013) (termed “Inverted Eclipse” therein) which are conjectured to arise from 
eclipses of a circumbinary disk. Other special variables of interest are also discussed. 

6 . We have found four previously unidentified eclipsing binaries, but it is not established whether 
they are members of the ONC or just line-of-sight associations. 
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A. Automatic and Subjective Selection Criteria 


Our original data had a varying range of photometric quality issues for a number of reasons; 
these quality issues were typically flagged in the reduced photometry catalog we retrieved from the 
WSA (see 12.2). These reasons include the bright, complex nebulosity in the ONC, the presence of 
many bright stars, the crowding towards the central Trapezium cluster, the fact that some stars are 
near chip edges, the variable weather from night-to-night, the fact that we were observing a huge 
number of stars, and the fact that some stars were near our saturation (brightness) or detection 
(faintness) limits for the exposure times we used. The process we used to select variables reflects 
these varying quality issues. Our criteria to bin stars into quality bins Q = 2, Q = 1, and Q = 0 is 
detailed below, and is primarily a function of how many bands (out of the three J, H, K bands) 
were deemed “pristine”, i.e., free of photometric error flags and within brightness and faintness 
limits that we considered safe. 


Q = 2 is reserved for stars with “pristine” data in all three bands. Stars were placed in the 
Q = 2 bin if they met all of these criteria: a) their PSF appeared starlike (i.e., not extended), 
b) they were observed between 50 and 125 times in each band . c) their mean J magnitude 


was between 11 and 17, H between 11 and 16, and K between 11 and 16, d) each band had no 
datapoints flagged for photometric quality issues, and e) they did not lie in the easternmost 3% of 
the observed field, which had anomalous data. 


Q = 1 is for stars with “pristine” data in one or two bands. Stars were placed in Q = 1 if they 
did not meet the Q = 2 criteria but a) their PSF appeared starlike (i.e., not extended), they did 
not lie in the easternmost 3% of the observed field, and one or two of their J, H, K bands met 
the following criteria: b) they were observed between 50 and 125 times in this band, c) their mean 
magnitude in this band was between 11 and 16 (for H and K) or between 11 and 17 (for J), d) no 
datapoints in this band had been flagged for photometric quality issues. 

Q = 0 is for stars with some suspicious or flagged data in all three bands. Stars were placed 
in Q = 0 if they did not meet Q = 2 oi Q = 1 but they were detected at least 50 times in at least 
one single band (i.e., at least 50 J or 50 77 or 50 K detections). 

Stars in Q = 2 and <5 = 1 were added to our final variable list if they met Stetson cuts, 
but no stars from the <5 = 0 bin were included in our final variable list without individual, by-eye 
scrutiny. Once stars are placed into the above quality bins, we select them as variable or nonvariable 
according to the following criteria. Note that Q = 2 and Q = 1 stars are selected automatically, 
while every Q = 0 star is classified by eye. Q = 2 stars are classified as variable if, in addition 
to meeting all of the Q = 2 criteria, their Stetson index exceeds 1.0. Q = 1 stars are classified as 
variable if their 5 > 1.0 and their one or two “pristine” data bands are variable independently of the 


^"^Notably, some stars lay on the overlap regions between detectors and were observed twice per night; sometimes 
these stars had erratic light curves due to a mismatch between the absolute flux calibrations of the two detectors 
(likely because they lay on an edge of one of these detectors). 
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rest. Q = 0 stars are classified as variable if they have S' > 1.0 and pass individual scrutiny where 
they are inspected for: a) bands co-vary in a consistent way through clear multiple-day trends, 
OR b) bands co-vary in a clearly periodic way. These two criteria are intended to exclude any 
variability that could feasibly be induced by instrumental or processing artifacts; as a side-effect, 
no short-term irregular variables were included via this criterion. 

To find periodic signals in Q = 0 stars, all Q = 0 stars with S > 1.0 were run through the 
period-finding search described in §3.2. 1[ Q = 0 stars that returned a strong period through this 
search were considered “subjective periodic candidates” and those that did not return a strong 
period were considered “subjective nonperiodic candidates”. 


The total pool of Q = 0 sources numbered 11,136. Of these, 2,273 had S > 1.0, though this 
Stetson value was influenced by corrupted data in a large number of cases. Our algorithm described 
in §3.2.1 detected periods in 231 of the Q = 0, 5 > 1 sources. After going through all of these 
2273 light curves by eye, we confirmed 153 non-periodic and 192 periodic Q = 0 stars for our final 
variable list. 


B. LOW-AMPLITUDE VARIABLE CANDIDATES 


While the focus of our present study is strong variables with S > 1.0, we examined possible 
variability in stars whose Stetson values lay between 0.55 and 1.0. The exact value of S which 
defines variability is difficult to assign, and depends on the coverage and noise properties of the 
data. For example. Carpenter et al. (2001) use S > 0.55 as the cutoff based on visual inspection of 
light curves along with concurrence from the value, which independently indicates that many 
of these stars are variable. Flaherty et al. (2013) used a similar method but selected S > 0.45 for 
their dataset. For the bulk of this study we focus on strongly variable stars with 5 > 1. However, 
some stars with Stetson index between 0.55 and 1.0 are certainly variable as well. We discuss this 
subsample as low-amplitude variable candidates, with the term “candidate” intended to convey 
that we may be seeing weak correlated instrumental or observational effects as well (such as were 
suspected by [Carpenter et n^|2001j). 


Figure shows the Stetson index S for all 2348 Q = 2 sources included in the initial study. 
Among these Q = 2 stars, 506 had a Stetson value within 0.55 < 5 < 1.0. We decided that, 
with only one exception, we could safely ignore all of these stars’ variability with no impact on our 
study’s objective. The single exception was a star whose very unusual anticorrelated variability 
between its J and K bands pushed its Stetson value down to S' = 0.97; we found its light curve at 
K band to have a period of 88.5 days after further analysis. We included this star as ONCvar 1226 
in our final variable list, and it is called out for discussion in §6.2| and §6.3[ 


The 506 low-amplitude variable candidate (1 < S < 0.55) light curves were period-searched 
in the same manner as described in (3.2.1 55 of them (10.8%) showed mathematically significant 
periods. This is roughly one-fourth the rate (42.1%) seen among Q = 2, S > 1 variables; these 
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disparate periodicity fractions may be an indication that many of the 1 < S' < 0.55 stars are 
not true variables. This is especially true since we expect most stars with such small variability 
amplitudes to be diskless stars, and diskless stars are noted in 13.3.2 to have a higher periodicity 
fraction than less-evolved, more stochastic higher-amplitude variables. 



